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SECTION  1 


INTRODUCTION 


Introduction 


Early  in  1972,  the  Advanced  Research  Projects  Agency  initiated 
a  several-year  program  to  develop  a  comprehensive  global  model  of  the 
ionosphere.  As  one  facet  of  the  overall  program,  Science  Applications, 

Inc.  (SAI)  emphasized  detailed  studies  of  the  chemistry  and  related  micro¬ 
scopic  processes  for  the  quiet,  mid-latitude  E-  and  F-regions  of  the 
ionosphere,  and  took  an  important  step  toward  cr  nsidering  high- latitude 
ionospheric  problems. 

SAI  first  modified  its  computer  code  (called  CHIEF)  for  chemistry 
of  ionization  in  the  E-  and  F-regions.  CHIEF  is  quite  elaborate,  carrying 
about  25  species  (including  excited  states)  and  several  hundred  reactions. 
This  code  was  originally  developed  for  studying  an  artificially-perturbed 
atmosphere.  The  modification  consisted  mainly  in  providing  the  code  with 
a  solar  source  of  UV,  EUV,  andx  rays  during  the  day  and  a  nighttime 
ionization  source.  As  a  step  toward  making  theoretical  predictions  with 
a  computer  program  about  the  chemical,  optical,  and  other  microscopic 
processes  in  a  natural  electron  aurora,  SAI  developed  an  electron  slowing- 
down  routine  —  with  account  of  discrete-energy  losses  instead  of  the  more 
usual  and  approximate  continuous  slowing-down  treatment  —  and  incor¬ 
porated  the  routine  into  the  large  ionospheric  chemistry  code,  now  called 
SOLAKEM.  This  code  was  tested  by  applying  it  to  the  experimental  con¬ 
ditions  appropriate  for  the  recently-reported  rocket  flight  of  Hays  and 
Sharp  [J.  Geophys.  Res.  78,  1153  (1973)]  who  have  made  the  most  accurate 
measurements  to  date  of  the  solar-produced  photoelectron-flux  spectra  in 
the  altitude  range  from  180  to  250  km.  SAI's  computed  photoelectron-flux 
spectra  are  in  very  good  agreement  with  the  Hays- Sharp  data.  The  details 
of  this  calculation  are  given  in  Section  2. 
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Additional,  lejs  extensive  studies,  also  to  be  published,  included 
(a)  determining,  for  a  12-species,  25-reaction,  E-  and  F-region  chemical 
kinetics  code,  the  uncertainties  in  species  densities  due  to  uncertainties 
in  the  rate  coefficients  for  the  reactant  species  (Section  3),  (b)  recognizing 
that  dissociative  recombination  of  molecular  ions  and  electrons  influence 
the  ionospheric  electron  temperature  to  a  larger  extent  than  some  of  the 
effects  usually  included  in  computing  the  electron  temperature  (Section  4), 
and  (c)  n  r.ing  that  the  temperature  dependence  for  the  dissociative  recom¬ 
bination  of  NO  customarily  used  in  E-  and  F-region  ionospheric  models 
is  incorrect  (Section  5). 
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CALCULATED  AND  OBSERVED 
PHOTOELECTRON-FLUX  SPECTRA  AT  DAWN 

i 

> 

> 


Calculated  and  Observed  Photoelectron-Flux  Spectra  at  Dawn 


B.  F.  Myers,  D.  A.  Hamlin,  and  M.  R.  Schoonover 
Science  Applications,  Inc.,  La  Jolla,  California  02037,  U.  S.  A. 

Abstract  —  Electron-flux  spectra  have  been  calculated  for  dawn  conditions 
at  altitudes  of  180,  200,  and  250  km  and  compared  with  recent  observations 
by  Hays  and  Sharp  (1973)  on  a  payload  launched  from  White  Sands,  New 
Mexico,  U.S.A.  on  February  8,  1971.  A  discrete  energy- loss  model  was 
used  in  the  calculations  and  good  agreement  is  found  between  the  calculated 
and  measured  spectra.  The  prominent  structural  features  in  the  observed 
spectra  are  reproduced  in  the  calculated  spectra  and  are  analyzed.  A  dis¬ 
crepancy  is  found  between  the  calculated  rate  of  electron  formation  by 
photoionization  and  collisional  ionization  and  the  loss  rate  deduced  from 
the  measured  electron  and  ion  densities.  No  definitive  explanation  of  this 
discrepancy  has  been  found,  either  in  terms  of  omitted  absorption  or 
transport  processes  or  in  terms  of  uncertainties  in  data  employed  in  the 
calculations. 


1.  INTRODUCTION 

Measurements  of  photoelectron- flux  spectra  have  recently  been  made  (Hays 
and  Sharp,  1973)  by  using  instrumentation  which  permits  meaningful 
measurements  in  the  range  of  electron  energies  from  below  one  electron 
volt  (eV)  to  several  tens  of  electron  volts.  These  flux  spectra  have  good 
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resolution  and  display  interesting  structural  features,  which  prompted  us 
to  test  a  method  of  calculating  the  electron-flux  spectra.  An  additional 
feature  of  the  experiment  was  the  simultaneous  measurements  of  the  charged 
particle  densities;  these  provide  further  tests  for  the  calculations. 

The  electron-flux  spectra  were  calculated  by  a  discrete  energy-loss 
procedure  with  an  energy  resolution  compatible  with  the  observations.  The 
procedure  is  part  of  a  larger  computer  code  which  is  used  to  describe  the 
chemical  and  energetic  transformations  in  the  E  and  F  regions  of  the 
ionosphere  over  short  time  intervals.  We  report  here  on  a  comparison  of 
the  calculated  and  observed  electron-flux  spectra  and  on  difficulties  in 
computing  energy  and  density  profiles  during  the  observations. 

2.  METHOD  OF  CALCULATION 

A  simple  method  was  used  to  treat  the  discrete  energy-loss 
processes  undergone  by  photoelectrons  formed  by  solar  radiation  incident 
on  atmospheric  species.  The  method  is  based  on  the  assumptions  that  a 
steady  state  lor  electrons,  sorted  into  a  set  of  energy  groups,  is  quickly 
established  after  photoabsorption  and  that  electrons  only  lose  energy  in 
each  interaction  process  until  a  specified  energy  limit  is  reached.  In 
the  present  application,  the  further  simplifications  of  local  energy 
deposition  and  isotropic  electron  fluxes  are  invoked.  The  total  set  of 
electrons  is  divided  into  a  ’’high- energy"  slowing-down  subset 


comprising  electrons  with  30  discrete  energies  between  1  and  103  eV  and 
a  "thermal”  subset  characterized  by  an  electron  temperature.  The 
relation  between  the  subsets  is  shown  in  Section  2.  2  to  be  in  accord  with 
a  recent  self- consistent  treatment  (Krinberg,  1973). 


2. 1  High-energy  electrons 

Application  of  the  steady-state  assumption  io  the  subset  of  high- 
energy  electrons  leads  to  the  following  expression  for  the  density  of  elec¬ 


trons  in  energy  group  j, 


[e.  (cm-3  bin"1)]  =  P  /L.  , 
J  J  J 


where  the  production  rate  is 


P.(cm"3  sec-1  bin  '‘)  =  EE  <p(i,t)  a(i,S)  b  (i,  j,S)  [S] 

J  o  :  P 


S  i 


photoabsorption 


+  E  i;.[ek]vkaxd  (k,S)bxd(k,j,S)[S] 
excitation,  dissociation,  and  dissoc.  excit. 


+  S  ^  [ek]  vk  o.(k,  S)  b.(k,  j,  S)  [S] 
ionization  and  dissoc.  ioniz. 

+  cj+i,j 


loss  to  thermal  electrons 


and  the  loss- rate  coefficient  is 

L.fsec'1)  =  Vj  {[1  -  bxd(i)j,S)]  oxda,S) 

+  [l-b1(j,i,S)]o.(j,S)}[S]  +  C.  .  ,  .  (3) 

1  '  J>J"A 

The  terms  on  the  right  hand  side  of  equation  (2)  refer,  respectively,  to 
the  processes  of  photoabsorption,  excitation  and  dissociation  (including 
dissociative  excitation),  ionization  (including  dissociative  ionization),  and 
loss  of  energy  to  thermal  electrons  by  the  high-energy  electrons;  the  latter 
three  processes  are  also  represented  by  terms  in  L.,  given  by  equation  (3). 
The  photoabsorption  term  gives  the  rate  of  forming  electrons  in  energy 
group  j  by  absorption  of  group-i  photons;  cp(i,t)  is  the  group-i  photon 
flux  at  time  t,  a(i,S)  is  the  photoionization  cross-section  for  species  S 
cbs orbing  group-i  photons,  and  b  (i,  j,S)  is  the  probability  that  a  group- j 

r 

electron  is  formed  when  species  S  absorbs  a  group-i  photon.  The  exci¬ 
tation  and  dissociation  term  gives  the  rate  of  forming  electrons  in 
group  j  by  group-k  electrons  impacting  on  S;  vfc  is  the  velocity  of  a 
group-k  electron,  °X£j(k>  S)  is  the  cross  section  for  excitation,  dissoci¬ 
ation,  or  dissociative  excitation  of  S  by  group-k  electrons,  and  b  (k,  j,S) 

All 

is  the  probability  of  forming  a  group-j  electron  when  a  group-k  electron 
interacts  with  S.  The  ionization  term  gives  the  rate  of  forming  group-j 
electrons  by  group-k  electrons  interacting  with  S;  c.(k,  S)  is  the  cross 
section  for  group-k  electrons  ionizing  S,  and  b.(k,j,S)  is  the  probability 


4 


r 


of  forming  group-j  electrons  by  group-k  electrons  interacting  with  S. 

f 

The  quantity  j  is  the  rate  at  which  electrons  leave  group  j+1  and 
enter  group  j  as  a  result  of  energy  loss  to  thermal  electrons.  Further 
t  consideration  of  C.+1^  as  well  as  of  the  probability  factors  bxd(j,  j,S) 

and  b.(j,  j,S)  of  equation  (3)  will  be  given  below. 

Ten  energy  groups  per  decade  of  electron  energy  were  established, 
as  shown  in  Table  1.  The  energy  groups  span  approximately  equal  inter¬ 
vals  on  a  logarithmic  scale;  AU.  is  the  width  of  group  j  and  U.  is  the 

J  J 

characteristic  energy. 

* 

The  resolution  of  the  energy- group  scale,  AU./U.,  is  slightly  above 

J  3 

20  per  cent  whereas  the  resolution  of  the  photoelectron  spectrometer  used 

k  in  the  observations  (Hays  and  Sharp,  1973)  being  modeled  was  10  per  cent. 

The  energy  resolution  in  the  model  could  be  increased  by  increasing  the 

number  of  groups,  but  comparison  of  calculated  and  observed  electron 
t 

fluxes,  as  shown  in  Section  4. 1,  indicates  no  need  for  such  change  in  the 
present  application. 

t  A  problem  arises  from  two  requirements  imposed  on  the  model. 

First,  the  model  must  accurately  describe  the  discrete  energy  losses  in 
relevant  physical  processes  and  second,  the  model  must  employ  an  arbi- 

t 

trary  set  of  energy  groups  and  related  characteristic  energies.  These 
requirements,  by  themselves,  are  not  compatible,  in  general,  since  the 
j  discrete  energy  loss  by  an  electron  in  a  particular  physical  process  is 

unlikely  to  have  a  value  equal  to  a  difference  in  characteristic  energies. 
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This  problem  is  overcome  by  ’’splitting"  electrons,  i.  e. ,  deriving  proba¬ 
bilities  that  the  outgoing  electron  in  the  physical  process  will  assume 
characteristic  energies  which  lie  on  either  side  of  the  true  value  of  the 
electron  energy.  The  probabilities  are  derived  on  the  basis  of  conserva¬ 
tion  of  energy  and  of  the  total  (free  plus  bound)  number  of  electrons. 

For  photoionization,  the  energy  of  the  outgoing  electron,  E  (i),  is 

G 

Ee(i)  =  Ep(i)  -  E(S')  ,  (4) 

where  Ep(i)  is  the  energy  of  the  gi  oup-i  photon  and  £(80  is  the  energy 
of  the  product  ion,  S'.  If 

U.  <  E0(i)  <  Uj+1  ,  j*l  (5) 

then  the  probabilities  for  the  electrons  are 


bp(i,j,S) 


b(i,  S,  S') 


(6) 


and 


bp(i,  j+1,  S)  =  Mi^SO  -  bp(i,j,S)  ,  (7) 

where  b(i,  S,S')  is  the  probability  that  product  S'  is  formed  by  target  S 
absorbing  a  group- i  photon.  This  probability  is  specified  by  experimental 


data  and  is  further  considered  in  Section  3. 1.  If,  instead  of  the  inequality 
of  equation  (5), 

Ee(i)  <  \J1  =  1. 1  eV  ,  (8) 

the  electron  is  added  to  the  thermal  subset  with  probability 
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b  (i,0,S)  =  b(i,S,SO 


(9) 


and  energy  E  (i)  *b(i,S,  S')  is  added  to  the  energy  of  the  thermal  subset. 

6 

For  an  excitation,  dissociation,  or  dissociative  excitation  event, 
the  procedure  for  deriving  probabilities,  t>xcj(k,  j,  S),  is  identical  to  the 
procedure  described  for  photoionization  events.  The  energy  of  the  out 
going  electron,  E  (xd),  is 

c 

E  ted)  =  U.  -  E(S')  (10) 

C  K 

where  U,  is  the  energy  of  the  impacting  electron  and  E(S')  is  the  energy 

XV 

of  the  product  S'.  Equations  (5)  through  (9)  apply  in  this  case  when  the 
quantities  E  (i),  b  (i,j,S),  b  (i,j+l,S),  and  bMjS')  are  replaced  by 

r  r  r 

Ee(xd),  bxd(k,j,S),  bxd(k,j+l,S),  and  bxd(k, S, S'),  respectively.  The 
quantity  bxd(k,S,S0  is  the  probability  for  forming  product  S'  by  a  group-k 
electron  interacting  with  target  S;  like  b(i,  S,  S'),  bxd(k,  S,  S')  also  is 
available  from  experimental  data.  In  calculations  of  b  (i,j,S)  and 

r 

bxd(k>  j,S),  the  numerical  accuracy  was  of  an  order  to  produce  errors  in 

-5 

energy  conservation  of  10  eV  or  less,  for  each  event. 

For  the  ionization  term  in  equation  (2),  a  single  ionization  event 
is  considered  in  which  there  are  two  outgoing  electrons;  the  lower- energy 
electron  is  named  the  "secondary”  and  the  other,  the  "degraded  primary" 
electron.  Both  electrons  are  described  by  differential  energy 
spectra  (Opal  et  al. ,  1971).  These  experimental  spectra  are  only 
approximated  by  the  procedure  employed  here  in  obtaining  values  of 
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bj(k,  j,S).  For  the  discrete  set  of  characteristic  energies,  XL,  and  energy- 
group  wiv'  «U.,  the  differential  probability,  i.  e. ,  probability  per  eV, 
for  forming  a  secondary  electron  with  specified  energy  is  chosen  to  satisfy 
the  relation 

imax  jmax 

§  Ps(i)'AUj  ■  Ps<W>  §  *Ui  s  Pa<W  W<W>  =  1  -  <«> 

where  the  differential  probability  for  the  secondary  electron,  P  (i  )  is 

sJmax 

assumed  to  be  constant  over  the  range  of  energies  assigned  to  the  secondary 
electrons.  ps(jmax)  is  restricted  to  certain  values  which  can  be  deduced 


from  equation  (11)  by  using  Table  1  for  values  of  W..  We  let  i  be 

]  Jmax 


determined  indirectly  by  choosing  the  largest  possible  value  of  P  (j  ) 

s  max 

that  is  less  than  the  differential  probability,  P,  for  the  secondary- electron 
energy,  Eg,  given  by 


when  evaluated  for  Eg  =  0.  The  functional  form  in  equation  (12)  is  that 
previously  used  by  Opal  et  al.  (1971)  to  describe  secondary- electron  spec¬ 
tra;  here,  we  slightly  generalized  by  replacing  the  ionization  potential  I 
by  I  +  X  where  X  is  the  electronic  excitation  energy  of  the  product  ion. 
The  quantity  E  is  a  shape  parameter  which  is  assumed  to  be  constant  for 
a  given  target,  as  it  varies  only  slowly  as  a  function  of  U.  ;  values  of  E 

K 

are  discussed  in  Section  3. 3. 
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The  differential  probabilities  defining  the  spectrum  of  degraded 
primary  electrons,  P'  and  P"  are  calculated  on  the  basis  of  energy 
conservation  and  unit,  total  probability;  two  differential  probabilities  for 
the  degraded  primaries  are  used  to  ensure  that  the  conservation  conditions 

can  be  satisfied  within  a  tolerably  small  error,  in  spite  of  the  fixed  values 
of  iU  .  Thus, 

J 


Jmax 

- 

r  n 

(U.-I-X)-P(j  ) 

K  s  Jmax 

2  AU.-U. 

j=c  J  J 

=  pd 
s 

£"Vuk 

+  pr 
d 
d 

\ 

(13) 


and 


- 

- 

X)au. 

k  k 

+  PH  * 

d'  d 

?AU* 

-  _ 

The  ranges  of  the  indices  k  and  l  are  selected  to  provide  positive  values 

of  Pd  and  and,  hopefully,  Pd  >  P';  the  inequality  is  in  accord  with 

observations  of  an  increasing  differential  probability  with  increasing  energy 

for  the  degraded  primary  electron  (Opal  et  al. ,  1971).  The  ranges  of  the 

indices  k  and  l  are  selected  by  an  algorithm  which  also  may  adjust  j 

Jmax 

to  satisfy  the  conditions  on  the  differential  probabilities  and  additionally 
select  the  maximum  value  of  P"  .  P'.  in  some  caseS)  the  inequality 

Pd  >  Pd  cannot  be  satisfied  and  in  these  cases  the  smallest  value  of 

Pd  ‘  Pd  is  chosen  trom  the  of  P',  P"  which  are  solution  sets  of 
equations  (13)  and  (14). 


Two  examples  of  the  use  of  the  method  for  computing  the  differ¬ 
ential  probability  spectra  in  a  single  ionization  event  are  showrn  in  Fig.  1; 
in  the  case  for  which  the  primary  energy  is  45  eV,  the  inequality  P''>  P' 
could  not  be  satisfied. 

The  relation  between  b.(k,  j,S)  of  equation  (2)  and  the  differential 
probabilities  Pg(jmax),  PJ,  and  P''  is 

b.(k,j,S)  =  P  AU  b(k,S,S')  (15) 

J  J 

where  P.,  denoting  Pg(jmax),  P^,  or  P",  is  the  differential  probability 
for  the  group  with  characteristic  energy  U.,  AU.  is  the  width  of  the  group, 
and  b(k,  S,S')  is  the  probability  for  forming  an  ion  S'  when  a  group-k 
electron  interacts  with  target  S;  the  probability  b(k,S,S')  is  available 
from  experimental  data  and  is  considered  in  Section  3. 3. 

The  probability  factors  in  equation  (3),  b  (j,  j,S)  and  b.(j,j,S), 
in  which  the  incident  and  an  outgoing  electron  are  assigned  to  the  same 
energy  group,  take  account  of  those  events  in  which  the  electron  energy 
loss  is  smaller  than  the  difference  -  Uj  These  factors  appear  in 

equation  (3)  as  a  result  of  factoring  the  density  [e.]. 

3 

The  loss-rate  coefficient  for  energy  loss  to  the  thermal  subset  of 
electrons  appearing  in  equations  (2)  and  (3)  is  calculated  with  the  approxi¬ 
mate  expression  (Schunk  and  Hays,  1971) 
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cj,]-l 


where 


7.73  x  10-6  [e  ] 
_ L  oJ 

6u.  (VuT) 

]  ] 


10  4 

2.  531  x  101U  U? 

_ _ I 


5.  387  x  1010  U. 


6U.  = 
3 


,  U.  <13.62  eV 
3 


,  U.>  13.62  eV 

J 


U.  -  U.  1  ,  j  a  2 
3  3-1  J 

U,  -  1  ,  j  =  l  , 


(16a) 


(16b) 


(  0.5  (Vu:  +yu7^)  ,  j  s  2 

<yu:>  = 

]  0.  5  (JU^  +  1)  ,  j  =  1  , 


(16c) 


and  [eQ]s  [e^]  is  the  density  of  electrons  in  the  thermal  subset.  The  form  of 
equation  (16)  is  valid  for  electron  temperatures  below  2500  °K.  For  the 
conditions  of  the  present  calculations,  the  energy  loss  per  collision  with 
the  thermal  subset  of  electrons  is  less  than  the  energy  difference  between 
groups;  this  is  the  basis  for  the  arguments  (j,j-l)  of  the  energy-loss 
coefficient  of  equation  (16). 


2.  2  Thermal  electrons 

The  thermal  subset  of  electrons  has  contributions  from  energy- 
degraded,  high-energy  electrons  as  a  result  of  the  loss  processes  included 


in  equation  (2),  as  well  as  from  photoabsorption  events.  The  rate  of 
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increase  in  the  thermal  electron  density  is  given  by  equation  (2)  with  j 
understood  to  represent  the  thermal  subset. 

The  total  electron  density  in  each  of  the  lower-energy  groups  of 
Table  1,  [e^.],  includes  a  significant  contribution  from  the  thermal  subset 
of  electrons  for  the  conditions  of  the  observations  being  modeled.  The 
fraction  of  the  thermal  electrons  in  group  j  is 


where  [e^]  is  the  density  of  thermal  electrons  in  group  j, 


(17) 


«W  =  erf(/x)--£r/F  e 

/? 


-x 


(18) 


erf(z)  is  the  error  function,  and 


yj  =  w./e  , 


(19) 


where  W.  is  the  upper-edge  energy  of  group  j  and  8  is  the  electron  tem¬ 
perature.  Equation  (18)  is  based  on  a  Maxwellian  distribution  of  energy, 


1(E) 


-  .2/1"  -E/e 


(20) 


J*  01 

The  total  electron  density  in  group  j  is  then 

[et,j]  =  tej]+  [eth,j]  •  (21) 

The  two  foregoing  paragraphs  describe  our  method  for  treating  the 
interaction  of  the  high-energy  electrons  with  the  ambient  electrons,  which 
is  in  good  agreement  with  a  rigorous,  self-consistent  treatment  (Krinberg, 
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r 


1973)  of  this  interaction.  To  illustrate  this  agreement  we  have  applied 
our  method  to  a  sample  problem  considered  by  Krinberg  (1973)  in  which 
the  photoelectrons  lose  their  energy  locally  by  only  Coulomb  collisions 
with  ambient  electrons;  the  photoelectron  production  spectrum  is 

q  ”  qo  eXp!'E  EV’  with  70  =  11)2  cm"3  sec'1  and  Eq  =  10  eV,  and  the 
electron  temperature  and  density  are  8  =  0.  25  eV  and  [eth]  =  5  X  lo5  cm'3, 

respectively.  For  this  sample  problem  the  production  rate  is 


P.(cm'3  sec'1bin'1)  =  q  (e  Wi_1  .  e"Wi)  *  r  r  l  , 

1  q°  6  )  +  Cj+l,j[ej+l]  •  W 


where 


w.  =  W./E  , 
J  J  o  * 

and  the  loss- rate  coefficient  is 


(22b) 


L.(sec_1)  --  c 
J  )  J,j-1  * 


For  the  slowing-down  electrons,  the  isotropic,  steradianal,  differential- 
energy  electron  flux  is 

<Mcm"2  sec"1  eV"1  sr"1)  =  m 

J  '  4«  ATT  »  (24 


4 it  AU. 
J 


where  v.  is  the  electron  velocity  corresponding  to  U..  The 
flux  for  a  Maxwellian  distribution  is 


corresponding 


v  f(E)  fe..  ] 
<ME)  =  - 


and  the  total  flux  is 


WSPPBWPWllfL,-  |j  .  .  J  I, .  I.  m  ,■ 4  *,..,4  ,...tp-,.i.  „  4  ju,ylt,i^,  t  j  n,.,  u  ^VM  _  ,,!>p*,  ,  .  ,.„j  w 


%i  *  *j  +  ■  (26) 

In  Fig.  2  the  triangular  points,  representing  the  total  flux,  were  computed 
from  equation  (26);  the  circular  points,  representing  only  the  slowing- 
down  flux,  were  computed  from  equation  (24).  The  curves  in  Fig.  2  are 
based  on  Krinberg's  (1973)  equation  (11),  modified  for  consistency  with 
lus  text;  the  solid  curve  represents  the  total  flux  and  the  dashed  curve, 
based  on  only  the  second  term  of  Krinberg's  (1973)  equation  (11),  neglects 
the  contribution  of  the  thermal  electrons  and  includes  only  the  slowing- 

down  flux.  The  agreement  between  the  two  methods  for  the  total  flux  is 
seen  to  be  quite  adequate. 

The  types  of  processes  undergone  by  the  thermal  subset  of  elec¬ 
trons  are  the  same  as  for  the  high-energy  subset,  except  for  energy 
limitations;  specific  processes  are  discussed  in  Section  3. 

2,3  Species  and  energy  densities 

As  a  result  of  the  processes  involving  the  photoelectrons,  species 
are  transformed  and  energy  is  redistributed.  Here,  these  changes  are 
distinguished  from  the  additional  changes  in  species  and  energy  densities 
involving  the  ambient  electrons,  chemical  reactions,  and  transport  of 
species  and  energy.  This  distinction  permits  some  understanding  of  the 
effect  of  the  formation  of  photoelectrons,  except  for  those  created  directly 
into  the  thermal  subset,  for  the  conditions  of  the  observations  (Hays  and 
Sharp,  1973).  The  rate  of  increase  in  the  total  electron  density, 
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[e]  =  teth]+  E  teji  ’ 

js  l  J 


(27) 


due  to  photoelectrons,  is 


ds] 

dt 


^(ij  t)  a(i,  S)  [S] 

S  ^  i  ^  ^ 

photoabsorption 

+  S  [ej]  v.  a.(j,S)  [S] 

ionization 


(28) 


The  rates  for  loss  of  species  S  and  gain  of  species  S',  due  to  photoelectrons, 


are 


d[S] 


dt  CT(i’S)  [S]  +  ^  [ej]  vj[oxd(j, S)  +  0.(j, S)]  [S] 


(29) 


d[S'] 


dt 


-  E  E  <p(i.  t)  c(i,  S)  b(i,  S,  S') 


S  i 


[S] 


+  ¥  j?i  [Si]  Vj[axd(i' S)  bxd(i- S-  S'>  +  VJ. S)  bi0’ s’ s')]  [S] 


(30) 


The  energy-density  changes  are  monitored  for  seven  modes: 
ionization,  dissociation,  eh  ctronic  excitation,  vibrational  excitation, 
radiation,  heavy-particle  translation  and  (molecular)  rotation,  and  elec¬ 
tron  translation.  The  heavy-particle  translational  and  rotational  modes 
are  assumed  to  be  in  equilibrium.  The  energy  transfers  into  the  ionization 
(I),  dissociation  (D),  and  electronic  excitation  (Xe)  modes  are  calculated  from 
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algebraic  expressions  such  as  ^  c  dS'/dt  where  c  represents  either  the 

S' 

ionization  potential  of  S',  the  dissociation  energy  per  atom  of  S',  or  the 


excitation  energy  of  S',  respectively,  and  dS'/dt  is  the  rate  of  increase 
in  S'  resulting  from  the  changes  being  distinguished.  Changes  in  the  re¬ 
maining  energy  modes  are  given  by  differential  equations  such  as 


dM 

dt 


M(i,S)  <p(i,t)  a(i,  S)  [S]  (1  -  6M  x  ) 


S  i 


M,  X, 


tej^l  vj[eMx^,S^  °xd^,S^  +  cMi^,S^  °iM(1  “  6M,Xv^S^ 

(31) 


-3, 


where  M(eV  cm  )  denotes  electron  translation  (E),  heavy-particle  trans¬ 
lation  (K),  radiation  (R),  or  vibrational  excitation  ( X y),  and  the  nega¬ 
tive  sign  is  used  only  for  M  =  E.  The  quantity  6^  ^  is  the  Kronecker 


delta.  In  equation  (31), 


=  E'MW.S'lW.S.SK'-VrJ 


r i  M 


M,X, 


(32a) 


Mx 


(i-s>  ■  5  Wi’s’s')bxd<i’s’s,) 


(32b) 


Wi-S>=  Vxv> 


(32c) 


where  cM(i,S,S'),  ^Mx(j, s, s'),  and  *Mi(j>s>s')  are  determined  by 
examining  the  energy  transformations  in  each  physical  process.  For 
M  =  R,  the  arguments  i  and  j  are  omitted. 
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2. 4  Procedure 


Calculations  of  electron-flux  spectra  and  short-time  energy 
and  density  profiles  are  made  with  our  E-  and  F- region  computer  code 
in  which  the  above-described  procedures  are  incorporated.  The  code  also 
computes  the  solar  flux  incident  on  the  volume  element  of  interest  by 
calculating,  with  appropriate  use  of  the  Chapman  function  (see,  e.  g. , 
Swider  u  v.  Gardner,  1969),  the  attenuation  of  85  photon  groups.  Inputs 
to  the  code  include  the  time  and  location  of  the  observations,  an  atmos¬ 
pheric  model,  and  species  and  energy  densities.  Optionally,  initial 
temperatures  may  be  used  in  place  of  the  electron,  heavy-particle,  and 
vibrational  energies,  from  which  the  initial  energies  are  calculated.  The 
changes  in  species  and  energy  densities,  and  then  temperatures,  are 
determined  by  integrating  differential  equations  which  include  the  terms 
of  equations  (28)  through  (31)  as  well  as  terms  accounting  for  many  other 
physical  and  chemical  processes.  Of  these  processes,  we  will  mention 
only  those  which  significantly  influence  the  results  discussed  herein. 

3.  PROCESSES  AND  DATA  SOURCES 

The  physical  and  chemical  processes  relevant  to  energy  loss  by 
photoelectrons  are  described  below.  Because  of  the  large  number  of  data 
employed  in  the  modeling  calculations,  only  data  sources  will  be  given, 


generally.  Lists  of  the  detailed  data  are  available  from  the  authors.  The 
cross-section  data  were  averaged  over  the  energy-group  widths  (Table  1) 
and  were  used  in  this  numerical  form  in  the  calculations. 

3. 1  Photoabsorption 

The  values  of  the  unattenuated  solar  fluxes  in  85  photon  groups  were 
taken  from  Hinteregger  (1970).  This  set  excludes  any  data  for  repeated 
spectral  intervals  of  Tables  I,  n,  and  in  and  the  first  four  entries  of 
Table  in  (Hinteregger,  1970).  A  slight  variation  in  these  data  was  also 
used,  as  discussed  below  in  Section  4. 1.  The  solar  photons  were  absorbed 
by  O,  Ng,  Og,  and  NO  in  their  ground  electronic  and  vibrational  states;  the 
range  of  wavelengths  for  the  solar  photons  included  is  30  to  1306  A. 

Absorption  by  0(  P)  was  limited  to  the  formation  of  the  five  atomic 
oxygen  ions  listed  in  Table  2.  The  cross  sections  were  chosen  from  the 
mean  of  theoretical  values  (Henry,  1967);  these  lie  between  the  experi¬ 
mental  data  (Cairns  and  Samson,  1965;  Comes  et  al. ,  1968).  Autoioniza¬ 
tion  (Henry,  1968)  was  not  included.  Production  distribution  probabilities, 
b(i,S,  S'),  were  obtained  from  Henry  (1967).  Instantaneous  radiation  by  the 
O  (  P)  and  O  (  P)  states  was  assumed;  see  Table  3. 

Absorption  by  Ng  resulted  in  the  products  listed  in  Table  2.  Cross 
sections  were  chosen  from  various  sources  (Huffman  et  al. ,  1963; 

Samson  and  Cairns,  1964;  Huffman,  1969;  Hudson  and  Carter,  1969; 
Huffman,  1972).  Recent  cross-section  data  (Carter,  1972)  of  generally 


higher  resolution  and  for  vibrationally- excited  ground- electronic-state  N0 

Ct 

(Cook  and  McNeal,  1972)  were  not  used  in  the  present  calculations  but  are 
considered  in  the  discussion  of  the  results  (see  Section  4.  2).  Use  of  inter¬ 
mediate  resolution  data  (Berkowitz  and  Chupka,  1969)  in  the  range  from 
700  to  785  A  would  not  significantly  change  the  adopted  cross-section 
values.  For  absorption  of  solar  photon  groups  representing  multiplets, 
the  cross  sections  for  lines  of  the  multiplets  were  averaged,  using  line 
strengths  (Wiese  et  al. ,  1966)  as  intensity  weighting  factors,  to  obtain  mul- 
tiplet  cross-sections.  In  determining  the  degree  of  ionization,  the  cross- 
section  data  references  given  above,  as  well  as  other  sources  (Samson  and 
Cairns,  1964;  Cook  and  Ogawa,  1965),  were  searched.  Where  the  degree 

of  ionization  was  less  than  unity,  the  formation  of  the  products  N(4S°)  and 
2 

N(  D°)  was  regarded  as  equally  probable.  At  wavelengths  above  the  ioniza¬ 
tion  threshold,  predissociation  was  considered  (Hudson  and  Carter,  1969). 

The  probabilities  for  ion  products,  b(i,  S,  S'),  were  obtained  from  Blake 
and  Carver  (1967)  for  wavelengths  of  580  A  and  greater;  at  smaller  wave¬ 
lengths,  the  probabilities  were  assumed  to  be  constant  at  the  values  for 
580  A.  For  the  ion  products,  only  the  ground  vibrational  level  was  assumed 
to  be  populated.  This  is  a  good  approximation  for  Ng+(X  2E  *)  (Blake  and 
Carver,  1967)  and  N„  (B  E  )  (Judge  and  Weissler,  1968),  but  for  N„+(A  2n  ), 
level  v'  =  1  may  be  significantly  populated.  Instantaneous  radiation  by 
Ng+(A  2nu)  and  Ng^B  2E^)  was  assumed;  see  Table  3. 
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The  11  products  of  absorption  by  Og(X  E^)  are  listed  in  Table  2. 
Cross-section  data  were  obtained  from  various  sources  (Blake  et  al. ,  1966; 
Cairns  and  Samson,  1965;  Huffman,  1969;  Huffman,  1972;  Matsunaga  and 
Watanabe,  1967;  Samson  and  Cairns,  1964;  Watanabe  et  al. ,  1953).  Cross 
sections  for  multiplet  photon- groups  were  obtained  as  for  Ng.  Photo¬ 
ionization  coefficients  were  obtained  from  Matsunaga  and  Watanabe  (1967). 
The  probabilities  b(i,S,  S')  at  wavelengths  between  the  ionization  threshold 
and  580  A  were  taken  from  Blake  and  Carver  (1967);  for  wavelengths 
smaller  than  580  A,  the  probabilities  were  assumed  to  be  constant  at  the 
values  for  580  A.  Further,  the  relation 

b  (i,02(X  1E^),  02+(a  4nu)j  =  4  b  (i,02(X  *Zg*),  02+(A  2IIu)j  (33) 

was  assumed;  however,  see  Carlson  and  Judge  (1971).  The  photodissoci¬ 
ation  products  are  not  well  known  (Matsunaga  and  Watanabe,  1967;  Beyer 
and  Welge,  1969;  Filseth  and  Welge,  1969).  Below  930  A,  estimates  of 
Matsunaga  and  Watanabe  (1967)  and  Beyer  and  Welge  (1969)  were  used;  at 
longer  wavelengths,  0(  P)  and  0(  D)  were  regarded  as  equally  probable 
dissociation  products.  As  in  the  case  of  Ng,  only  the  ground  vibrational 
level  of  ions  was  assumed  to  be  populated;  this  is  not  a  good  assumption 


for  oxygen,  but  apparently  it  is  not  serious  for  the  results  of  the  present 

+  4  -  +2 

calculation.  Instantaneous  radiation  by  Og  (b  E^  ),  Og  (A  fl^),  and 
+  2  - 

On  (  E  )  was  assumed;  see  Table  3. 

2  g 
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The  products  of  absorption  by  NO  were  restricted  to  the  list 
of  six  in  Table  2.  Cross-section  data  were  taken  from  several  sources 
(Cook  and  Ching,  1965;  Sullivan  and  Holland,  1966;  Watanabe,  1958; 
Watanabe  et  al. ,  1967).  Resolution  of  multiplets  was  not  employed  for 
these  data  as  for  Ng  and  Og.  Ionization  coefficients  (Cook  and  Ching, 

1965;  Watanabe  et  al. ,  1967)  were  arbitrarily  set  to  unity  below  580  A. 

The  product  distribution  was  estimated  cn  the  basis  of  phase-space 
considerations. 

3.  2  Collisional  excitation  and  dissociation 

3 

The  four  products  of  excitation  of  0(  P)  are  given  in  Table  2; 
cross-section  data  for  each  product  state  were  taken  from  Henry  et  al. 

(1969)  and  Stone  and  Zipf  (1971).  Prompt  radiation  by  0(3p  ^P)  and  0(3s  ^S° 
was  assumed  (see  Table  3). 

See  Table  2  for  the  10  products  considered  in  excitation  of 

N0(X  +).  For  excitation  to  the  A  +,  B  ^11  ,  C  ,  and  E  **£  + 

2  g  u  g  u  g 

states,  data  of  Borst  (1972),  Borst  et  al.  (1972),  and  Jobe  et  al.  (1967) 

were  used  but  with  the  recognition  of  potential  difficulties  (Shemansky  and 

Broadfoot,  1971a,  1971b).  For  energies  larger  than  40  eV,  these  data 

-3 

were  extrapolated  by  using  an  E  -dependence  of  the  cross  section,  as 

determined  theoretically  for  triplet  states  (Chung  and  Lin,  1972).  The 

cascades  B  -»  A  +,  C  -»  B  ,  and  E  +  ->  A  +  (however, 
g  u’  u  g’  g  u  v 

see  Brinkman  and  Trajmar  (1970)),  were  assumed  to  occur  instantaneously. 
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The  excitation  cross-sections  to  the  a  n  state  for  E  <  40  eV  (Borst, 

b 

1972)  and  for  2000  *  E(eV)  *  100  (Holland,  1969)  are  highly  compatible; 

cascade  to  the  a^n  state  is  small  (Borst,  1972;  Ajello,  1970).  The  radi- 

6 

ative  decay  of  the  a  1I1  state  was  assumed  to  occur  instantaneously;  see 

s 

Table  3.  For  excitation  to  the  c'  *E  *  state,  we  used  cross-section  values 
equal  to  one  half  of  those  of  the  theoretical  calculations  (Chung  and  Lin, 
1972)  as  a  compromise  between  theory  and  experiment  (Aarts  and  De  Heer, 
1971),  and  assumed  that  the  dipole- allowed  transition  c'  *11  +  ->  X  ^E* 
dominates.  For  excitation  to  the  a*  *E  +  state;  we  used  the  data  of 

B 

Brinkmann  and  Trajmar  (1970);  radiative  decay  to  the  ground  electronic 

state  of  N_  was  assumed  for  the  conditions  of  the  observations  being 
z 

modeled,  in  spite  of  the  facts  that  the  transition  is  only  quadrupole-allowed 

(Dressier  and  Lutz,  1967)  and  other  decay  modes  need  to  be  considered 

(Brinkmann  and  Trajmar,  1970).  Excitations  '.o  the  dipole- allowed  states 

b'  *E  +  and  b  *11  were  included  by  using  theoretical  cross-section  curves 
g  u 

(Chung  and  Lin,  1972).  For  the  b  *11  state,  1/10  of  the  computed  cross- 

section  was  used  since  this  state  is  derived  from  a  mixed  configuration  and 

the  cross  section  has  been  computed  for  only  one  of  two  configurations; 

this  is  by  analogy  with  the  case  of  the  b'  *E  +  state  where  taking  mixing 

into  account  reduces  the  pure-configuration  cross-sections  by  a  factor  of 

10.  The  W  Au  state,  as  well  as  others,  was  neglected.  However,  for 

the  W  3A  state,  at  least,  recent  observations  (Chutjian  et  al. ,  1973) 
u 
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yield  substantially  larger  cross-sections  than  predicted.  For  vibrational 

excitation  of  N0(X  *E  +),  we  used  the  data  of  Chen  (1966)  for  electron  en- 
2  g 

ergies  up  to  3.  5  eV  and  other  sources  for  larger  energies  (Comer  and 
Read,  1971;  Pavlovic  et  al. ,  1972).  Interpolation  was  required  in  the 
energy  range  between  3.  5  and  18  eV. 

The  dissociative  excitation  of  Ng(X  E^),  as  treated  here,  leads 
to  the  formation  of  the  eight  states  of  atomic  nitrogen  listed  in  Table  2, 
with  the  product- partner  N(  S°),  followed  by  instantaneous  decay  as  given 
in  Table  3.  Cross-section  data  for  each  state  were  taken  from  several 
sources  (Sheridan  et  al. ,  1961;  Aarts  and  De  Heer,  1971;  Mumma  and  Zipf, 
1971)  as  compiled  by  Kieffer  (1972),  but  some  extrapolation  of  the  data  was 
required,  particularly  at  the  higher  electron  energies.  The  treatment  of 
dissociative  excitation  is  also  somewhat  hampered  by  lack  of  knowledge  of 
intermediate  states  in  the  excitation  processes. 

3  - 

The  nine  products  of  excitation  and  dissociation  of  Og(X  Eg)  are 
given  in  Table  2.  For  excitation  to  the  a  A^  and  b  Eg  states,  the  data 
(T  raj  mar  et  al.,  1971)  were  extended  to  energies  above  45  eV  by  assuming 
an  E" ^-dependence  of  the  cross  section.  The  c  *£u -state  excitation  data 
(Trajmar  et  al. ,  1972)  are  limited  and  were  extended  by  interpolation  and  ex¬ 
trapolation.  The  c  *£  “  state  of  Og,  which  is  not  carried  as  a  species  in  the 
E-  and  F-region  code,  is  assumed  to  dissociate  immediately  after  formation 
and  yield  two  0(3P)  atoms.  For  the  B  3eJ  state,  excitation  is  similar  to  that 


of  the  c  “  state.  The  data  (Trajmar  et  al. ,  1972)  are  limited  and 
estimates  (Peterson  et  al. ,  1969)  have  been  made;  the  guessed  cross- 
section  curve  includes  data  on  excitations  at  9.  97  and  10.  29  eV.  Dissoci¬ 
ation  to  0(3P)  +  0(1D)  was  assumed  to  be  the  dominant  result  of  the  processes 
represented  by  the  guessed  curve.  Data  on  vibrational  excitation  is  limited 
(Linder  and  Schmidt,  1971)  to  electron  energies  below  2  eV,  although 
direct  excitation  to  high  vibrational  levels  has  been  observed  (Trajmar  et 
al. ,  1971).  In  the  present  calculation,  vibrational  excitation  was  included 
only  for  electron  energies  less  than  2  eV. 

For  dissociative  excitation  of  02(X  )  to  the  3s  S°  state  of  atomic 

oxygen  (Ajello,  1970;  Lawrence,  1970;  Mumma  and  Zipf,  1971;  Aarts  and 

De  Heer,  1971)  and  to  the  3s5S°  state  (Wells  et  al.  ,  1071;  Wells  and  Zipf, 

3 

1972),  the  product  partner  was  taken  to  be  0(  P).  Tue  excited  states  of 
atomic  oxygen  were  assumed  to  radiate  instantaneously;  see  Table  3.  Other 
excited  states  produced  by  dissociative  excitation  of  02(X  ‘  Eg  )  (Judge, 

1972)  were  neglected. 

For  excitation  of  0(3P)  and  excitation  and  dissociation  of  N2(X  Eg  ) 
and02(X  3£  “),  the  probabilities  bxd(k,S,S')  were  determined  from  the 
set  of  cross  sections  for  all  processes  included  in  this  work. 

3.3  Collisional  ionization 

The  total  ionization  cross-section  data  for  O^P)  at  electron  ener¬ 
gies  greater  than  40  eV  (Fite  and  Brackmann,  1959;  Rothe  et  al. ,  1962) 


and  below  40  eV  (Boksenberg,  1961)  as  compiled  by  Kieffer  (1969),  were 
combined.  The  probabilities  b.(j,  S,  S')  were  taken  from  data  of  Seaton 
(1959)  for  electron  energies  below  250  eV;  above  250  eV,  the  probabilities 
were  taken  to  be  constant  at  the  values  corresponding  to  250  eV. 

For  ionization  of  Ng(X  ),  the  total  cross-section  data  used  were 
those  of  Rapp  and  Englander- Golden  (1965).  By  using  these  data  together 
with  other  sources  in  the  case  of  formation  of  the  A  2IIu  state  (Shemansky 
and  Broadfoot,  1971b;  Holland  and  Maier,  1973),  of  the  B  2£u+  state  (Borst 
and  Zipf,  1970),  and  of  the  dissociative  ionization  (Rapp  et  al. ,  1965;  also 
see  Kieffer  and  Van  Brunt,  1967),  we  estimated  the  cross  section  for  pro- 
duction  of  the  X  state.  The  data  on  formation  of  B  +  were  multiplied 
by  a  factor  of  1. 12  to  account  for  bands  not  included  in  measurements.  The 
radiative  decay  processes  of  the  A  and  B  2Eu+  states  were  assumed  to 
be  instantaneous;  see  Table  3.  The  C  2£u+  and  D  2n  states  were  neglected. 
Also,  the  process  of  dissociative  ionization  was  greatly  simplified  by  as¬ 
suming  that  the  products  of  this  process  are  the  neutral  and  ionized  nitrogen 
atoms  in  their  ground  electronic  states;  note  that  a  reasonably  good  as¬ 
sumption  is  that  no  radiative  states  of  N  II  are  formed  in  the  dissociative 
ionization  procesr  (Aarts  and  De  Heer,  1971). 

The  six  states  included  in  the  ionization  of  09(X  ~)  are  listed 

2  g 

in  Table  2.  Total  ionization  cross-section  data  were  obtained  from  Rapp  and 
Englander- Golden  (1965);  data  for  other  states  are  not  extensive  in  the  oxygen 


4  - 

system.  For  excitation  to  the  b  L  state,  total  cross-sections  were 

S 

calculated  from  the  available  data  (Borst  and  Zipf,  1970).  Adjustments  to 

o 

the  data  (Skubenich,  1968)  for  excitation  to  the  A  state  were  made  by 

(a)  multiplying  the  data  by  the  ratio  of  the  result  of  Borst  and  Zipf  (1970) 

4 

to  that  of  Skubenich  (1968)  for  excitation  to  the  b  ZJ  state  and,  (b)  at 
energies  above  the  limit  of  the  data,  140  eV,  using  a  constant  fraction, 

n 

0.35,  of  calculated  cross-sections  (Peterson  et  al. ,  1969)  for  the  A  n 
state.  The  fraction  0.35  gave  agreement,  on  the  average,  between  calcu¬ 
lated  and  adjusted  data  below  140  eV.  This  fraction  was  also  used  to  adjust 

4 

calculated  values  (Peterson  et  al. ,  1969)  for  excitation  of  the  a  Ilu  state. 

2 

The  cross  sections  for  excitation  to  the  II  state  were  obtained  as  the 

g 

residual  values  of  the  total  cross-section  data  (Rapp  and  Englander- Golden, 

1965)  after  accounting  for  other  excitations  and  dissociative  ionization 

(Rapp  et  al. ,  1965).  The  relative  values  of  the  II  -state  excitation  cross- 

g 

sections  obtained  by  this  procedure  are  somewhat  larger  than  calculations 
predict  (Peterson  et  al. ,  1969);  clearly,  additional  data  are  needed,  but 
in  the  present  calculations  this  lack  is  not  as  critical  as  it  is  in  general, 
due  to  the  smaller  contributions  from  processes  involving  molecular 
oxygen  at  the  altitudes  of  interest.  The  recognized  prompt  radiative  transi¬ 
tions  following  ionization  of  molecular  oxygen  are  given  in  Table  3.  The 
treatment  of  dissociative  ionization  in  molecular  oxygen  was  analogous 
to  that  for  molecular  nitrogen. 
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For  ionization  processes,  the  probabilities  b.(j,S,S0  were  derived 
in  the  same  way  as  for  excitation  and  dissociation  processes.  The  values 
of  E  for  ionization  of  Ng  and  02>  as  used  in  equation  (12),  were  taken  from 
Opal  et  al.  (1971);  E(N2)  =  12. 7,  E(02)  =  17. 4.  For  O  and  the  dissociative 
ionization  of  N2  and  C>2,  E  was  estimated  by  assuming  that  E  s:  I  +  X. 

3-  4  Energy  quantities 

Energy  quantities  are  required  in  equations  (4),  (10),  (31),  and  (32). 
These  were  obtained  by  using  established  values  of  the  term  energies 
(Table  2)  and,  for  dissociative  excitation  or  ionization,  by  using  the  values 
of  the  kinetic  energy  of  the  heavy-particle  fragments.  For  dissociative 
ionization  of  N2(X  L+),  roughly  3  eV  per  atom,  when  this  energy  is  avail¬ 
able,  is  attributed  to  the  atomic  fragments  (Kieffer  and  Van  Brunt,  1967; 
also,  see  Stockdale  and  Deleanu,  1973).  The  same  energy  per  atom  was 

used  for  dissociative  excitation  of  ftyX  h  +)  and  dissociative  ionization 
3  ® 

of  02(X  Lg  )  in  the  absence  of  experimental  data.  In  dissociative  excita- 
3 

ti°n  o£  °2(X  Eg  )>  rouShly  0.  5  eV  per  atom,  when  available,  is  attributed 
to  the  product  atoms  (Borst  and  Zipf,  1971;  also,  see  Freund,  1971; 
Stockdale  and  Deleanu,  1973). 

The  energy  quantities  for  the  diatomic  products  are  for  the  ground 
vibrational  level  of  the  electronic  states  listed;  thus,  except  for  the  ground 
electromc  states  of  N2  and  C>2;  no  account  was  taken  of  the  vibrational 
excitation  which  occurs  in  the  excitation  and  ionization  of  diatomic  species. 
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3.  5  Thermal  electrons 


The  thermal-subset  electrons  undergo  some  of  the  above-described 
processes  and  others  in  addition.  These  interactions  are  important  for 
their  contributions  to  the  electron  densities  in  the  low-j  groups  (see  Table  1) 
and  to  the  total  electron  density  and  energy- density  changes.  The  inter¬ 
actions  of  thermal  electrons  will  be  briefly  described. 

For  dissociative  recombination  of  NO+(X  *E+),  experimental  data 
(Weller  and  Biondi,  1968)  have  been  combined  with  theory  (Bardsley,  1970) 
by  using  an  inverse  square-root  dependence  on  the  electron  temperature,  0. 
The  effect  ohvery  recent  measurements  (Walls  and  Dunn,  1974)  on  the  dis¬ 
sociative  recombination  of  NO+,  not  available  when  the  present  calculations 
were  made,  is  discussed  in  Section  4.  2.  Combined  data  (Kasner  and 
Biondi,  1968;  Smith  and  Goodall,  1968;  Mehr  and  Biondi,  1969;  Mahdavi 

et  al. ,  1971;  also,  see  Cunningham  and  Hobson,  1972)  were  used  with  a 

-0  7  +2 

temperature  dependence  of  0  *  for  dissociative  recombination  of  Og  (X  n 

and  Og+(a  Ilu).  For  Ng+(X  ^g+)>  data  of  Mehr  and  Biondi  (1969)  were  ad¬ 
justed  to  pass  through  the  300°K-value  of  Kasner  (1967). 

Thermal  electrons  undergo  elastic  collisions  with  ions  (Banks,  1966) 
and  neutrals  (Banks,  1966).  (For  electron-neutral  elastic  collisions,  the 
numerical  values  derived  (Valerio,  private  communication)  and  used  here, 
from  equation  (20)  together  with  equations  (25),  (27),  and  (29)  of  Banks 
(1966),  are  smaller  by  a  factor  of  1.  50  than  those  given  in  Table  3  of  Banks 


g' 
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(1966)  for  Ng,  Og,  and  O. )  We  also  included  rotational  excitation  and 

de-excitation  of  Ng(X  *E*)  and  Og(X  ^E  ~)  (Dalgarno,  1969);  vibrational 

1  +  3  - 

excitation  and  de-excitation  of  N0(X  E  )  (Chen,  1966),  00(X  E  )  (Linder 

^  S  4  g 

o 

and  Schmidt,  1971),  and  of  NO(X  n)  (Spence  and  Schulz,  1971);  excitation 
and  de-excitation  among  the  ground- configuration  terms  of  neutral  and 
ionized  atomic  oxygen  and  nitrogen  (Seaton,  1955;  Saraph  et  al. ,  1966; 
Czyzak  et  al. ,  1967;  Moiseiwitsch  and  Smith,  1968);  collisional- radiative 
recombination  with  atomic  ions  of  oxygen  and  nitrogen  (Pates  et  al. ,  1962; 
Bates  and  Kingston,  1964);  and  cooling  of  thermal  electrons  by  excitation 
of  the  atomic-oxygen  ground-term  levels.  For  this  fine- structure  cooling 
we  used  the  results  of  Breig  and  Lin  (1966),  following  the  procedure  of 
Dalgarno  and  Degges  (1968),  with  account  of  radiative  transfer  (Craig  and 
Gille,  1969);  however,  doubts  about  the  treatment  of  this  mechanism  of 
electron  energy  loss  remain  (Takayanagi  and  Itikawa,  1970). 

Transport  processes  involving  (essentially)  thermal  electrons  are 
discussed  in  Section  4.  2. 

3.  6  Observational  conditions 

The  observations  being  modeled  were  obtained  (Hays  and  Sharp,  1973) 
on  the  flight  of  NASA  rocket  13.51  launched  February  8,  1971  at  White  Sands, 
New  Mexico  (L  32.°413  N,  X  106.°322  W)  at  6:56:18  LT.  From  more  extensive 
data  (Sharp,  private  communication),  we  computed  zenith  angles,  including 
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an  approximate  correction  for  the  rocket  trajectory.  The  calculations  of 
electron-flux  spectra,  described  in  Section  4. 1,  were  made  at  180- ,  200- , 
and  250- km  altitude,  for  which  the  calculated  zenith  angles  are  90?  17, 

90?  12,  and  89?  95,  respectively.  Small  errors  may  exist  in  these  values 
as  a  result  of  incomplete  data,  but  these  possible  errors  have  negligible 
effect  on  the  results  and  discussion  (Section  4). 

The  zenith  angles  were  used  in  the  E-  and  F-region  code  to  compute 
the  solar  flux  incident  at  the  altitudes  of  interest  in  each  of  the  85  solar 
photon  groups.  Attenuation  was  attributed  to  the  species  0(  P),  Ng(X  £g), 
and  0„(X  3E  ");  the  densities  of  these  species  as  well  as  the  temperature 
of  the  neutral  species  were  obtained  from  two  atmospheric  models. 

(a)  Jacchia  (1971),  with  an  exospheric  temperature  of  1000° K  (Hays  and 
Sharp,  1973)  and  (b)  essentially  that  used  by  Hays  and  Sharp  (1973)  but 
extended  to  include  Ng.  At  120- km  altitude,  we  assumed  [02]  =  4  x  1010 
cm'3,  [Nj  =  4  x  1011  cm"3  (CIRA,  1965),  and  tho  ratio  [0]/[02]  =  3.  5  as 
in  Hays  and  Sharp  (1973).  At  higher  altitudes,  we  calculated  densities  and 
temperatures  according  to  equations  (4)  and  (2)  of  Walker  (1965),  with 
the  exospheric  temperature  To  =  1000°K  and  c  =  0.01617  (Jacchia,  1965; 
Fig.  1,  curve  2b). 

The  initial  densities  of  the  ions  and  electrons  required  by  our 
E-  and  F-region  code  were  taken  from  the  flight  measurements  (Hays  and 
Sharp,  1973).  An  important  minor  species  is  NO;  since  no  measurements 


exist  for  the  flight  conditions,  we  used  model  calculations  of  the  density 
(St  rob  el,  1972).  Densities  of  other  species  included  in  the  chemistry  code 
are  unimportant  for  the  present  calculations.  Initial  values  of  energy 
densities  were  calculated  from  electron  temperatures  (Hays  and  Sharp, 
1973),  neutral-particle  temperatures  (discussed  above),  particle  densities, 
and  the  assumption  of  equal  neutral-particle  and  vibrational  temperatures. 

4.  RESULTS  AND  DISCUSSION 

The  principal  results  of  the  calculations  reported  here  are  electron- 
flux  spectra  at  altitudes  of  180,  200,  and  250  km  for  the  conditions  of  NASA 
rocket  flight  13.  51  on  February  8,  1971.  These  results  are  compared  with 
measurements  of  the  electron-flux  spectra.  The  effect  on  the  calculations 
of  restrictions  imposed  by  the  flight  conditions  on  the  rates  of  change  of 
charged- particle  densities  is  discussed.  In  addition,  photoabsorption  and 
photoelectron  interactions  are  described  in  terms  of  formation  of  species 
and  partition  of  energy  for  this  flight. 

4. 1  Electron-flux  spectra 

The  calculator  electron-flux  spectra  at  three  altitudes  are  compared 
with  the  measured  spectra  in  Fig.  3.  The  general  result  is  one  of  good, 
apparent  agreement  between  the  calculated  and  measured  spectra;  how¬ 
ever,  to  make  such  comparisons  for  spectra  at  altitudes  much  above 
200  km  is  questionable.  Recent  calculations  (Mayr  et  al. ,  1973)  for 


energy-deposition  profiles  show  that  at  250  km  about  1/5  or  more  of  the 
energy  is  being  deposited  over  an  altitude  range  greater  than  the  scale 
height  of  atomic  oxygen.  (This  result  is  only  indicative  of  the  problem,  as 
the  photoelectron  spectrum  and  possibly  other  conditions  of  the  calculations 
(Mayr  et  al. ,  1973)  are  not  precisely  applicable  here.)  Thus,  the  circum¬ 
stances  of  the  agreement  at  250  km  of  the  calculated  and  measured  fluxes 
apparently  differ  from  those  at  the  lower  altitudes;  nevertheless,  this  pos¬ 
sibility  will  be  put  aside  in  some  qualitative  discussions  below.  The  com¬ 
plete  list  of  calculated  flux  values  is  given  in  Table  4. 

The  effect  of  uncertainties  in  the  input  data  on  the  electron-flux 
spectra  has  been  examined  in  only  a  very  limited  way  by  varying  the  atmos¬ 
pheric  composition  and  the  solar-flux  values.  The  use  of  the  atmospheric 
model  of  Jacchia  (1971),  on  which  the  calculated  electron  fluxes  in  Fig.  3 
are  based,  yields  electron-flux  values  which  are  smaller  by  less  than 
20  per  cent  than  those  values  obtained  by  using  our  extended  version  of  the 
model  of  Hays  and  Sharp  (1973).  Two  solar-flux  changes  have  been  ex¬ 
amined.  If  all  the  solar-flux  values  are  changed  by  a  factor  up  to  2.  5, 
then  the  electron-flux  values  are  changed  by  roughly  the  same  factor  in  the 
energy  range  where  thermal  electron  contributions  are  negligible.  If  the 
solar-flux  values  recommended  by  Donnelly  and  Pope  (1973)  are  adapted  to 
the  photon  group  structure  of  the  present  calculations,  the  electron-flux 
values  differ  from  those  displayed  in  Fig.  3  by  an  average  of  less  than 
20  per  cent  at  180  km  and  by  smaller  amounts  at  higher  altitudes.  The 
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differences  result  predominantly  from  the  lower  flux  values  of  Donnelly 
and  Pope  (1973)  in  the  100  to  200  A  range. 

The  uncertainty  in  the  measured  values  of  the  electron-flux  spectra 
at  250  km  is  indicated,  for  random  errors,  on  Fig.  3.  At  lower  altitudes 
the  counting  rates  are  smaller  and  the  errors  are  larger  (Sharp,  private 
communication);  the  counting  rate  is  roughly  a  factor  of  three  lower  at 
180  km  than  at  250  km. 

There  are  several  prominent  features  of  the  measured  electron-flux 
spectra  which  we  discuss  in  terms  of  the  calculated  spectra.  The  well- 
known  minimum  between  2  and  3  eV,  resulting  from  the  electron-induced 
vibrational  excitation  of  Ng,  is  evident  for  altitudes  of  180  and  200  km 
(Fig.  3).  As  the  altitude  increases,  this  minimum  vanishes  since  the  nitro¬ 
gen  density  decreases  and  the  production  of  electrons  with  energies  in  this 
range  increases  as  the  electron  density  increases.  The  contribution  of 
thermal  electrons  to  the  calculated  flux  in  the  2-  to  3-eV  range  of  electron 
energy  is  negligible  at  altitudes  of  180  and  200  km,  but  at  250  km  only  8  per 
cent  of  the  flux  is  from  the  high-energy  subset  of  electrons.  Two  discrep¬ 
ancies  may  be  noted  at  180  and  200  km.  First,  the  calculated  minimum  is 
much  deeper  than  the  measured  one  and,  second,  the  minimum  occurs  in 
the  energy  range  from  2.  0  to  2.  5  eV  in  the  calculations  but  in  the  range 
from  2.  5  to  3.  2  eV  in  the  observations.  The  first  discrepancy  would  not  be 
rignificantly  reduced  by  accounting  for  de-excitation  of  vibrationally- excited 
nitrogen  by  high-energy  electrons  in  energy  groups  below  2  eV  (in  the 
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thermal  subset  of  electrons,  both  excitation  and  de- excitation  are  included); 
inclusion  of  vibrational  de- excitation  would  increase  the  fluxes  by  less  than 
10  per  cent  and  20  per  cent  for  the  energy  ranges  from  2. 0  to  2.  5  eV  and 
from  2.  5  to  3.  2  eV,  respectively.  Furthermore,  significant  error  in  the 
total  cross-section  for  vibrational  excitation  of  Ng  seems  unlikely,  based 
on  recent  experiments  (Spence  et  al. ,  1972).  By  contrast,  the  cross  sec¬ 
tions  for  excitation  of  the  ground- configuration  terms  of  atomic  oxygen  are 
uncertain;  these  processes  are  important  in  the  electron  energy  range  in 
question,  and  only  theoretical  values  are  available.  However,  doubling  or 
halving  the  cross  section  does  not  change  the  flux  value  in  the  range  from 
3. 2  to  4.  0  eV  at  180  km  (Fig.  3),  and  doubling  the  absolute  cross-sections 
while  maintaining  relative  values,  reduces  the  depth  of  the  minimum  (at 
3.  2  eV  on  Fig.  3)  by  only  10  per  cent  while  increasing  the  difference  in 
flux  values  at  4.0  eV.  Thus,  unless  the  shape  of  the  theoretical  cross- 
section  curve  (Henry  et  al. ,  1C39)  is  in  serious  error,  the  uncertainty  in 
cross-section  values  for  excitation  of  the  ground- configuration  terms  of 
atomic  oxygen  will  apparently  not  account  for  the  first  discrepancy. 

The  second  discrepancy,  in  the  energy- location  of  the  minimum, 
is,  from  the  viewpoint  of  the  calculations,  a  result  of  the  larger  cross- 
section  for  vibrationally  exciting  Ng  in  the  range  from  2.  0  to  2.  5  eV  than 
in  the  range  from  2.  5  to  3.  2  eV;  a  better  energy  resolution  in  the  calcula¬ 
tions  might  shift  the  minimum  but  this  has  not  been  determined.  A  greater 
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electron  temperature  would  remove  the  discrepancy  but  the  minimum 
value  required,  1352  K,  falls  outside  the  apparent  range  of  scatter  of  the 
measured  temperatures  (Hays  and  Nagy,  1973)  at  180  km. 

A  possibility  for  altering  both  discrepancies  simultaneously  is  the 
inclusion  of  a  photoelectron  flux  from  the  region  conjugate  to  the  observa¬ 
tion  region.  The  former  region  is  sunlit  in  the  present  case.  If  the  con¬ 
jugate  photoelectron  flux  is  as  large  as  that  derivable  from  the  measure¬ 
ments  of  Maier  and  Rao  (1972),  then  the  depth  of  the  calculated  minima  for 
180  and  200  km  would  be  reduced  by  a  factor  of  two  roughly  and  the  minima 
would  be  shifted  into  the  2.  5  to  3.  2  eV  energy  range.  For  larger  energies, 
the  qualitative  aspects  of  the  electron-flux  spectra  would  be  unaltered  and 
flux  values  would  be  increased  by  less  than  50  per  cent  on  the  average.  By 
contrast,  if  the  conjugate  photoelectron  flux  is  as  small  as  that  calculated 
by  Nagy  and  Banks  (1970;  see  also  Nagy  et  al. ,  1973),  then  the  calculated 
electron-flux  spectra  at  180  and  200  km  (Fig.  3)  would  be  essentially 
unaffected. 

A  plateau  occurs  between  20  and  30  eV  in  the  electron-flux  spectra 
at  200  and  250  km  for  both  the  calculations  and  observations,  as  shown  in 
Fig.  3.  This  plateau  is  a  result  of  the  absorption  of  the  He  II  solar  line 
at  303.  8  A  by  O  and  N2;  Fig.  4  shows  a  comparison  between  the  electron 
fluxes  with  and  without  the  contribution  from  photoabsorption  of  the 

O 

303.  8  A-He  II  line.  The  increased  fluxes  in  the  energy  group  from  20  to 
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25  eV,  shown  in  Fig.  4,  are  due  principally  to  the  formation  of  electrons 

in  photoionization  processes  in  which  the  excited  electronic  states  of  the 

3  1  + 

ionized  target  are  produced;  for  absorbers  0(  P)  and  N9(X  E  ),  the 
principal  heavy- particle  products  are  0+(^D°),  0+(^P°),  N~(A  ),  and 

Z  11 

_j_  2  + 

Ng  (B  Lg  ).  This  phenomenon  occurs  also  at  180  km,  but  there  the 
increased  flux  from  absorption  of  the  303.  8  A-He  II  line  amounts  to  a 
change  of  only  1. 4  per  cent  at  23  eV  and  does  not  form  a  distinct  plateau. 
With  increasing  altitude,  photoelectron  production  by  absorption  of  the 
303.  8  A-He  H  line  increases  significantly  and  results  at  250  km  in  a  rela¬ 
tively  larger  electron  flux  in  the  23  eV  energy  bin  than  in  the  adjacent 
bins.  Over  the  range  of  energies  from  14  to  36  eV,  the  loss-rate  coeffi¬ 
cient  L.  in  equation  (1)  has  relative  values  at  180,  200,  and  250  km  which 
1 

change  by  not  more  than  15  per  cent;  the  contribution  of  O  dominates  in 
the  loss- rate  coefficient  at  these  altitudes  and  this  fact  is  responsible  for 
the  small  change  in  the  relative  values  of  the  loss- rate  coefficient  in  spite 
of  the  density  changes  in  O  and  Ng. 

The  observations  yield  only  hints  of  a  ledge  above  35  eV  in  the 
electron-flux  spectra;  the  calculations  predict  a  plateau  in  the  energy  range 
above  35  eV  as  shown  in  Fig.  5.  This  ledge  results  from  the  photoabsorp- 

O 

tion  of  solar  radiation  in  the  wavelength  range  from  about  153  to  205  A. 


This  fact  is  demonstrated  in  Fig.  5  for  three  altitudes  by  comparison  of 
the  electron  flux  calculated  with  and  without  the  contribution  of  solar 
radiation  in  the  stated  wavelength  interval.  According  to  the  data  of 
Hinteregger  (1970),  which  are  used  here,  this  wavelength  interval  is 
divided  at  176  A  into  two  subintervals  for  which  only  the  integral  flux  is 
presented  without  further  assignment  of  the  radiation  to  specific  sources 
and  wavelengths;  however,  further  assignments  can  be  made  (Donnelly 
and  Pope,  1973).  Nevertheless,  the  principal  contribution  to  the  ledge 
is  from  photoabsorption  processes  in  0(3P)  and  N9(X  *£  +)  for  which  the 
product  ions  are  in  electronically-excited  states:  0+(3D°),  0+(3P°), 

O  (2s2p  P),  Ng  (A  nu),  and  Ng  (B  ^u+).  The  decreasing  slope  of  the 

ledge  with  increasing  altitude  is  related  to  the  photoabsorption  of  greater 
wavelength  radiation  which  produces  electrons  in  the  energy  group  with 
U16  =  36  eV;  this  contribution  increases  more  rapidly  with  increasing 
altitude  than  does  the  contribution  of  photoelectrons  to  the  energy  group 

with  =  45  eV  from  photoabsorption  of  radiation  in  the  wavelength  range 
from  153  to  205  A. 

The  energy  resolution  in  the  present  calculations  seems  adequate 
for  modeling  the  observations,  as  indicated  in  Fig.  3.  Nevertheless, 
with  increased  energy  resolution,  further  structure  in  the  electron-flux 
spectra  could  be  expected  to  appear  and  is  observed  in  very  recent 
observations  (Mukai  and  Hirao,  1973).  To  bring  the  present  calculations 


into  perspective,  the  spectrum  of  the  differential  photoelectron  production 
rate  for  the  energy  interval  from  20  to  30  eV  is  shown  in  Fig.  6;  this 
spectrum  is  based  on  the  photoabsorption  data  described  in  Section  3. 1  but 
is  averaged  over  0. 2-eV  intervals.  On  the  expectation  that  the  structure 
of  this  spectrum  would  be  reflected  in  an  electron-flux  calculation  in  which 
the  energy  resolution  was  of  the  order  of  tenths  of  an  electron  volt,  a 
qualitative  comparison  of  the  spectrum  with  the  observations  (Mukai  and 
Hirao,  19/3)  is  encouraging  for  further  calculations,  in  spite  of  the  differ¬ 
ent  conditions  of  the  calculations  here  and  the  observations  of  Mukai  and 

Hirao  (1973).  A  similar  qualitative  comparison  is  made  by  Mukai  and 
Hirao  (1973). 

4-2  Electron  production  and  loss 

While  the  calculated  and  observed  electron-flux  spectra  are  in 
reasonable  agreement,  there  is  an  underlying  inconsistency  between  the 
production  of  electrons  by  photoionization  and  collisional  ionization  and  the 
loss  of  electrons  by  dissociative  recombination.  Since  the  rocket  flight 
(Hays  and  Sharp,  1973)  was  at  dawn,  the  net  rate  of  change  in  electron 
density  is  expected  io  be  positive.  However,  the  loss  rate  exceeds  the  pro¬ 
duction  rate  by  factors  of  15,  10,  and  4.  8  at  altitudes  of  180,  200,  and 
250  km,  respectively.  See  Table  5.  The  precise  values  of  these  factors 
depend  on  our  choice  of  data.  For  data  of  one  type,  such  as  the  unattenu¬ 
ated  solar  fluxes,  an  alternate  choice  of  values  could  reduce  these  factors 
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but  at  the  expense  of  degrading  the  agreement  between  the  calculated  and 
observed  electron-flux  spectra.  For  data  of  another  type,  alternate  choices 
would  predominantly  change  only  the  stated  factors.  Contributions,  includ¬ 
ing  those  associated  with  the  latter  type  of  data,  which  may  be  important 
for  the  inconsistency  will  now  be  examined.  In  the  following,  "discrepancy" 
will  refer  to  the  positive  difference  between  the  electron  loss  and  produc¬ 
tion  rates;  while  several  of  the  possible  contributions  can  significantly 
reduce  the  discrepancy,  no  definitive  explanation  for  the  discrepancy  has 
been  found. 

A  possible  source  of  error  in  calculating  the  electron  loss  rate  is 
the  uncertainty  in  the  temperature  dependence  of  the  dissociative  recom¬ 
bination  rate  coefficient.  This  is  most  likely  to  be  so  for  the  dissociative 
recombination  of  NO+  where  the  only  previously  available  guide  to  the 
correct  temperature  dependence  for  ionospheric  conditions  was  theory 

(Bardsley,  1970).  On  this  basis,  in  modeling  the  Hays  and  Sharp  (1973) 

-1/2 

observations,  we  assumed  the  rate  coefficient  varies  as  8  .  Very 

recently,  however,  ion-storage  techniques  (Walls  and  Dunn,  1974)  have 

been  used  to  obtain  cross  sections  for  the  dissociative  recombination  of 

several  molecular  ions  as  a  function  of  electron  energy.  For  NO+,  these 

measurements,  when  combined  with  low-temperature  experiments  (Weller 

and  Biondi,  1968),  indicate  a  stronger  dependence  of  the  rate  coefficient 

-1/2 

on  the  electron  temperature  than  0  .  As  an  approximate  upper  limit 
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to  the  enhanced  dependence,  consider  a  dependence  of  0  .  Use  of  this 

dependence  reduces  the  electron  loss  rate  by  dissociative  recombination 
by  roughly  40  per  cent.  But  further  theoretical  calculations  (Michels, 

1973)  predict  a  dependence  of  0“°*  56  for  temperatures  below  2000°K,  so 
that  consistency  in  knowledge  about  the  temperature  dependence  of  the 
dissociative  recombination  rate  coefficient  is  apparently  not  yet  achieved. 

Another  possible  error,  uncovered  by  the  work  of  Gutcheck  and 
Zipf  (1973),  is  that  previous  experiments  (Bardsley  and  Biondi,  1970)  about 
the  dissociative  recombination  of  NO+,  02+,  andN2+  may  have  involved 
vibrationally- excited  molecular  ions,  and  this  fact  may  have  resulted  in 
derived  rate  coefficients  which  are  in  error.  For  NO+,  02  ,  and  Ng  ,  the 
rate  coefficients  most  probably  decrease  as  the  vibrational  energy  level 
of  the  reacting  molecular  ion  increases  (O’Malley,  1969);  for  dissociative 
recombination  of  N2+,  this  is  apparently  so  only  for  v  ^  2  (except  in  the 
temperature  range  between  680  and  3900°K  where  the  rate  coefficient  for 
v  =  2  exceeds  that  for  v  =  0)  (Michels,  1972).  Thus,  if  the  experiments 
(Bardsley  and  Biondi,  1970)  did  involve  vibrationally- excited  molecular 
ions,  correction  of  the  analyses  of  these  experiments  would  likely  yield 
larger  rate  coefficients  and  consequently  increase  the  discrepancies  calcu¬ 
lated  in  the  modeling  of  the  Hays  and  Sharp  (1973)  observations.  However 
note  that  in  the  case  of  dissociative  recombination  of  C>2  ,  shock-tube 
experiments  (Cunningham  and  Hobson,  1972)  and  ion-storage  experiments 


(Walls  and  Dunn,  1974)  are  essentially  in  agreement  with  afterglow  experi¬ 
ments  (Weller  and  Biondi,  1968)  and  superficial  considerations  indicate 
that  conditions  in  the  three  experiments  are  not  likely  to  be  similar  in 
regard  to  vibrational  excitation  of  the  molecular  ions.  Further  work  is 
needed  to  establish  the  occurrence  of  errors  in  the  rate  coefficients  as  a 
result  of  the  possibility  undercovered  by  Gutcheck  and  Zipf  (1973),  but,  at 
present,  this  possibility  cannot  be  used  to  account  for  the  discrepancies  in 
electron  production  and  loss  rates  under  the  conditions  of  the  observations 
(Hays  and  Sharp,  1973). 

In  the  above  argument,  the  tacit  assumption  was  made  that  excited 
vibrational  states  of  molecular  ions  are  not  important  under  quiet  iono¬ 
spheric  conditions;  this  is  probably  true  for  NO+  (Myers,  1973)  but  has 
not  been  investigated  for  O*  and  N2+. 

A  further  error  source  in  computing  the  electron  loss  rate  is 
associated  with  uncertainties  in  the  measured  densities  of  electrons  and 
molecular  ions;  the  errors  here  were  not  explicitly  given  (Hays  and  Sharp, 
1973),  so  further  discussion  will  be  left  to  others. 

Contributing  to  errors  in  calculation  of  the  rate  of  change  in  elec¬ 
tron  density  would  be  omission  of  transport  processes  for  electrons  and 
ions.  Consider  the  transport  of  charged  species  by  ambipolar  diffusion, 
by  interaction  with  neutral  winds  and  electric  fields,  and  by  atmospheric 
expansion.  Ideally,  the  contribution  of  diffusion  (Schunk  and  Walker,  1973) 


could  be  calculated  from  the  experimental  data;  however,  since  this  con¬ 
tribution  depends  principally  on  the  values  of  the  second  derivatives  of 
density  and  temperature  with  respect  to  altitude,  the  calculation  is  not 
practical  when  based  on  experimental  data,  as  a  result  of  the  uncertainties 
in  measured  temperatures  (Hays  and  Nagy,  1973)  and  presumably  densi¬ 
ties.  However,  under  the  restrictions  of  a  monotonic  increase  in  electron 
density  and  electron  and  ion  temperatures  with  increasing  altitude  and  of 
positive  values  for  the  second  derivatives  of  the  density  and  temperatures, 
one  can  demonstrate  that  diffusion  is  unlikely  to  account  for  as  much  as  a 
third  of  the  discrepancy  and,  more  probably,  a  considerably  smaller 
fraction.  In  this  estimate,  the  second  derivatives  of  the  density  and  tem¬ 
peratures  had  fractional  changes  with  altitude  similar  to  those  of  the  cal¬ 
culated  discrepancies  over  the  range  from  180  to  250  km. 

To  determine  the  contribution  of  neutral  winds  to  the  rate  of  elec¬ 
tron  density  change  is  even  more  difficult  than  that  of  diffusion.  There 
were  no  measurements  of  winds  associated  with  the  rocket  flight  (Hays 
and  Sharp,  1973)  and  appropriate  calculations  of  the  winds  have  not  been 
made.  However,  models  of  neutral  winds  (Kohl  and  King,  1967;  Sterling 
et  al. ,  1969;  Strobel  and  McElroy,  1970)  indicate  that  during  the  night,  a 
wind  component  results  in  flow  from  north  to  south  and  thereby  causes  an 
upward  drift  of  ionization  along  geomagnetic  field  lines,  whereas  during  the 
day,  the  wind  direction  and  ionization  drift  are  reversed.  The  models  also 


indicate  that  during  the  time  of  the  rocket  flight,  the  north-south  compo¬ 
nent  of  the  wind  had  just  or  was  about  to  change  direction,  with  the  conse¬ 
quence  that  velocities  were  of  the  order  of  tens  of  meters  per  second  rather 
than  of  the  order  of  hundreds  of  meters  per  second  which  nay  occur  at 
other  times.  Furthermore,  if  neutral  winds  are  to  contribute  to  removing 
'screpancies  of  the  present  calculations,  a  positive  value  of  the  rate 
of  electron  density  change  must  obtain.  On  the  basis  of  these  considera¬ 
tions,  an  estimate  indicates  that  the  interaction  of  neutral  winds  with 
Charged  species  accounts  for  less  than  20  per  cent  of  the  discrepancies. 

The  contribution  of  electric  fields  is  also  probably  unimportant,  as 
calculations  have  demonstrated  (Raster,  1971)  that  neutral  winds  are  more 
important  than  electric  fields  in  their  effect  on  the  behavior  of  the  F  region. 

Finally,  the  expansion  of  the  atmosphere  (Sterling  et  al. ,  1969), 
which  has  a  nearly  maximum  rate  during  the  time  of  the  rocket  flight,  re¬ 
sults  in  an  electron  loss  rate  which  is  less  than  1.  5  per  cent  of  the 
discrepancies. 

Accounting  for  electrons  from  the  region  conjugate  to  the  observa¬ 
tion  region  (see  Section  4. 1)  could  increase  the  calculated  electron  produc¬ 
tion  rate  and  thus  decrease  the  discrepancies.  But,  conjugate  electron 
fluxes  based  on  recent  measurements  at  the  altitudes  of  interest  (Maier 
and  Rao,  1972)  or  on  a  combination  of  measurements  and  calculations  (Nagy 

and  Banks,  1970;  Nagy  et  al. ,  1973)  would  reduce  the  discrepancies  by  less 
than  6  per  cent. 


43 


Errors  in  calculating  the  photoelectron  production  rate  can  arise 
from  certain  aspects  of  the  present  treatment  of  the  photoionization  processes 
and  from  omission  of  processes.  (Below,  ’’photoelectron  production”  will 
refer  to  formation  of  electrons  by  only  photoionization. )  In  certain  spec¬ 
tral  intervals,  average  absorption  cross-s  actions  were  used  where  there 
is  significant  variation  of  the  cross  section  with  wavelength.  The  use  of 
average  cross-sections  leads  to  errors  in  calculating  the  incident  solar 
flux  at  the  altitudes  of  interest  and  to  an  additional  error  in  calculating  the 
amount  of  photoabsorption.  At  photon  energies  hv  less  than  the  thresholds 
hi^Ng)  and  hv^Og)  for  ionization  of  Ng  and  Og,  photoionization  of  Og  and 
NO  will  be  treated  incorrectly;  at  hv  >  [hv.(N0),  hv.(09)],  the  amount  of 
photoionization  of  Ng  and  O ^  will  be  in  error.  The  effect  of  using  average 
absorption  cross-sections  was  examined  for  Ng  at  hv(eV)  <  ^(Ng)  +  1  and 
for  Og  at  hv(eV)  <  hv^Og)  +  1  by  recomputing  the  photoionization  of  Ng, 

Og,  and  NO  with  more  accurate  cross-section  values  or  selected  cross- 
section  values  which  maximized  the  amount  of  photoionization  in  the  spectral 
interval  of  averaging.  For  hv  <  [hv^Ng),  hv.COg)],  the  increase  in  photo¬ 
ionization  of  Og  and  NO  was  less  than  25  per  cent  of  the  presently  computed 
rate  of  photoelectron  production  at  180  km  under  the  conditions  of  the 
observations  (Hays  and  Sharp,  1973)  and  consequently  less  than  1  per  cent 
of  the  discrepancy  between  electron  loss  and  production  rates.  For 
hv>  [hv^Ng),  hv^Og)],  the  change  in  the  amount  of  photoioniz^.tion  between 
180  and  250  km  was  negligible  with  respect  to  the  discrepancies. 
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The  photoelectron  production  rate  may  be  increased  by  including 
vibrationally-excited  absorbers  with  smaller  ionization  thresholds  than  for 
unexcited  absorbers.  We  estimated  the  effect  of  including  vibrationally- 
excited  N2  in  the  photoabsorption  processes  by  (a)  using  absorption  data 

for  vibrationally-excited  N2  (Cook  and  McNeal,  1972),  (b)  assuming  a 

vibrational  temperature  of  2500%  which  is  larger  than  expected  (Breig 
et  al. ,  1973),  and  (c)  omitting  attenuation  of  the  solar  radiation.  The 
increase  in  the  photoelectron  production  rate  by  N2  absorbing  radiation 
in  the  range  from  796  to  842  A  was  less  than  12  Per  cent  of  the  rate  without 
considering  vibrationally-excited  N2  and  was  negligible  in  comparison  with 
the  discrepancies.  For  vibrationally-excited  02,  no  photoionization  data 
are  available  in  the  range  from  1027  to  1080  A;  however,  for  an  assumed 
vibrational  temperature  of  2500%  the  contributions  of  vibrationally- 
excited  02  to  the  photoelectron  production  rate  would  be  negligible  in  com¬ 
parison  with  the  discrepancies,  provided  the  absorption  cross-sections 
are  less  than  about  10'17  cm2.  But  for  computing  the  photoelectron  pro¬ 
duction  rate,  these  contributions  could  be  non-negligible  provided  the 
photoionization  cross-section  of  vibrationally-excited  02  is  effectively 
5  x  10-18  cm2  or  larger  for  the  C  II  and  O  VI  multisets  in  the  vicinity  of 
1032  and  1037  A,  and  the  solar  flux  values  of  Donnelly  and  Pope  (1973) 
are  essentially  correct. 
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The  effect  of  rotationally-excited  Ng  and  O 2  is  also  negligible  with 
regard  to  the  electron  production  and  loss  discrepancies. 

Certain  absorptions  by  atomic  oxygen  (Huffman  and  Larrabee,  1968) 
were  neglected  in  our  calculations.  Of  these  absorptions,  only  that  of  the 
H  Lyman-beta  line  by  0(  P2)  is  non- negligible.  On  the  basis  of  the  ratio 
of  the  Lyman- alpha  flux  to  the  Lyman-beta  flux  observed  at  approximately 
185  km  (Detwiler  et  al. ,  1961),  the  attenuation  of  the  Lyman-beta  flux  in 
our  electron-flux  calculations  is  in  error  to  the  extent  of  overestimating 
the  photoelectron  production  rate  by  about  4  per  cent;  this  is  small  com¬ 
pared  with  the  discrepancies  at  180  km.  At  higher  altitudes  this  error  in 
attenuation  is  even  less  important. 

The  uncertainty  in  [NO]  could  significantly  affect  the  discrepancy 
at  180  km,  listed  in  Table  5,  since  the  photoionization  of  NO  accounts  for 
about  25  per  cent  of  the  photoelectrons  formed  at  180  km.  The  density  of 
NO  was  based  only  on  calculations  (Strobel,  1972). 

Under  other  conditions,  Swartz  and  Nisbet  (1973)  have  found  a 
discrepancy  similar  to  the  one  discussed  above  and  Roble  and  Dickinson 
(1973)  have  found  a  related  one.  These  investigators  suggest  that  under¬ 
estimation  of  solar  fluxes  might  account  for  the  discrepancies;  this  possi¬ 
bility  appears  less  tenable  in  the  present  calculations  for  the  reason  stated 
at  the  beginning  of  this  section.  Increases  in  the  solar  EUV  flux  values  by 
factors  of  two  to  three  as  proposed  (Swartz  and  Nisbet,  1973;  Roble  and 
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Dickinson,  1973)  would  leave  the  present  calculation  with  a  still  significant 
discrepancy  between  electron  production  and  loss. 

4*3  Species  and  energy  partitioning 

The  species  formation  and  the  transformation  of  energy  associated 
with  photoabsorption  and  the  subsequent  interactions  of  the  high-energy 
electrons  have  been  compiled  for  the  conditions  of  the  observations.  As 
stated  in  Section  2.3,  this  compilation  excludes  changes  involving  ambient 
electrons,  chemical  reactions,  and  transport  of  energy  and  species.  In 

Table  6,  the  relative  formation  rates  for  species  are  listed  when  these 
rates  equal  or  exceed  0. 01  of  the  total;  the  total,  absolute  rates  are  273, 

378,  and  535  cm'3  sec'1  for  altitudes  of  180,  200,  and  250  km,  respectively. 
Several  of  these  values  are  noteworthy.  The  large  fractional  contribution 
for  0(  D)  arises  from  electron  collisional  excitation.  The  fractional 
contribution  of  0(3s  3P)  seems  rather  large  in  view  of  its  excitation  energy 
but  this  is  due  to  a  combination  of  two  factors;  a  sharply  peaked  cross- 
section-versus-energy  curve  for  formation  of  0(3s  3P)  by  electron  impact 
in  the  energy  range  between  12  and  32  eV  (Stone  and  Zipf,  1971)  with  a 
large  maximum  cross-section  and  a  fractional  increase  in  electron  flux 
with  increasing  altitude  which  is  largest  in  the  same  energy  range.  For 
formation  of  vibrationally  excited  nitrogen,  N2(X  \Y ,  the  maximum 
as  a  function  of  altitude  in  the  relative  rates  follows  from  the  maximum  in 
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the  product  of  densities  for  electrons  and  nitrogen  molecules,  the  electron 
density  increasing  and  the  nitrogen  density  decreasing  with  altitude.  The 
formation  of  Ng(A  2S+)  is  predominantly  by  absorption  of  radiation  with 
energy  less  than  the  ionization  threshold  for  nitrogen  as  treated  here.  A 
more  careful  treatment  of  this  process  would  undoubtedly  reduce  the  listed 
relative  rates  of  formation  and  increase  those  of  atomic  nitrogen  species. 

The  energy  partitioning  is  shown  in  Table  7  in  terms  of  the  relative 
rates  of  change  in  the  seven  energy  modes.  See  Section  2. 3.  For  each 
altitude,  the  first  column  shows  the  partitioning  immediately  following 
photoabsorption  (PA)  and  the  second  column  shows  the  partitioning  when 
the  steady-state  in  the  high-energy  electron  subset  has  been  established 
(SS),  i.  e. ,  after  this  subset  has  excited,  dissociated  and  ionized  the  target 

species. 


5.  CONCLUSIONS 

Electron- flux  spectra,  calculated  for  180- ,  200- ,  and  2 50- km 
altitude,  are  in  reasonable  agreement  with  electron-flux  spectra  measured 
for  dawn  conditions  in  midlatitude  (Hays  and  Sharp,  1973).  The  agree¬ 
ment  indicates  that  the  processes  and  data  employed  are  sufficient  for 
calculations  with  spectral  energy  resolution  of  approximately  20  per  cent. 
Certain  prominent  spectral  features  are  shown  to  be  associated  with  photo¬ 
absorption  processes  in  which  electronically- excited  ions  are  formed; 


greater  spectral  resolution  in  both  calculations  and  observations  is  needed 
to  probe  more  deeply  into  these  spectral  features. 

The  loss  rate  of  electrons  by  dissociative  recombination,  based 
on  measured  electron  and  ion  densities,  exceeds  the  calculated  rate  of 
forming  electrons  by  photoionization  and  collisional  ionization.  This  dis¬ 
crepancy  decreases  with  increasing  altitude.  No  definitive  explanation 
for  the  discrepancy  was  found,  although  a  number  of  processes  associated 
with  photoabsorption,  particle  transport,  and  data  errors  were  examined. 
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Table  1.  Energy  groups  used  in  treating  electron  energy  loss  processes 


Group  No.  ,* 

j 

VeV 

t 

AU.,  eV 

t  ^ 

w  =  £  AU 

J  i=0  1 

0 

0.  5* 

1.0 

1.0 

1 

1.1 

0.25 

1.25 

2 

1.4 

0.35 

1.6 

3 

1.8 

0.4 

2.0 

4 

2.3 

0.5 

2.  5 

5 

2.9 

0.7 

3.2 

6 

3.6 

0.8 

4.0 

7 

4.5 

1.0 

5.0 

8 

5.7 

1.4 

6.4 

9 

7.2 

1.6 

8.0 

10 

9.0 

2.0 

10.0 

11-20 

10  uHo 

10  AUj-10 

10  Wj-io 

21-30 

100  Uj-20 

100  AUi-20 

100  Wj— 20 

The  high-energy  subset  of  electrons  are  those  electrons  in  group  j=l 


lTj  is  the  characteristic  energy  of  group  j;  AU.  is  the  width  of  group  j: 
w.  is  the  upper- edge  energy  of  group  j.  ■* 

vVhen  it  is  necessary  to  adopt  probabilities  for  an  electron  to  enter 
both  the  thermal  subset  and  group  1,  a  characteristic  energy  of 
0.  5  eV  is  assigned  to  the  lower- energy  electron,  which  energy  is  then 
added  to  the  total  thermal-electron  energy. 
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Table  2.  Species  and  Energy  Quantities 


O^P) 

a,xd,  I 

o(4D) 

xd 

a,xd 

0(JS) 

xd 

a 

0(3s  5S°) 

xd 

0(3s  V) 

xd 

a,xd 

0(3p  5P) 

a 

0(3p  3P) 

xd 

a 

o*;4s°) 

a,  1 

1 

o+(2d°  ) 

a,i 

o+(2p°) 

a,l 

0+(2s2p4  4P) 
0+(2s2p4  2P) 

0,(xV,  v>0) 

°2(a  JV 

°2<«  At> 

°f  P 

°2*<X,V 
°2>  "«> 
°2<\"u> 

°2  <S  V 

°2  <  V 


a 

a 


xd 
xd 
xd 
xd 
xd 
a,i 
a.l 
a,  I 
a,i 
a 


2.  5G1 
1.96 
4.17 
9. 14 
9.52 
10. 74 
10.99 

13.62 
16.  93 

18.62 
28.49 
39.98 

0. 98 
1. 63 
4.05 
6.12 
12.06 
16.10 
16.81 
18.2 
20.0 


N(3s  4P)  xd 

N(3s  2P)  *  xd 

N(2p4  4P)  xd 

N(3p4D°)  xd 

N(3pV)  xd 

N(3s'  2D)  xd 

N(4s  2P)  xd 

N(3d  2D)  xd 

N2(x1E^,v>0)  xd 

N2(A!Su)  a'xd 

N2(D3ng)  xd 

N2(a  !ng)  xd 

N2(C3n  )  Xd 

N2(E  V ■)  xd 

N2(a'  V)  xd 

N2(h  \)  Xd 

N2(c'  \+)  xd 

N2(b'  \+)  xd 

N+(3P)  i 

N2+(XV>  M 

N2+(A2nu)  a,l 

N2+(B20  a,l 

NO+(X  V) 


a 


4.881 
2.37 
3.56 

10.33 
10. 67 
10.  93 
11.70 
11.78 

12.34 
12.86 
13.02 

6. 17 
7.35 
8.  55 
11.02 
11.87 
12.  28 

12.57 
12.  93 
13.21 
14.53 

15.58 
16.  58 
18.75 
10. 20  ’ 


•a  .  photoabsorptlon;  ad  »  colll.lonat  citation  and  dissociation;  i  .  coltlsional  ionization 

»J2,rSlSMe.r,W  u,oho„o”Crn  0<’P)  or  Nl"S>  ““P*  •«  those  Uttar; 
Energy  dillorcnccs  among  levels  o  o  3p‘  t,  Jl'  .Z  15  I'™”"1"1  hy  02(x  3;  ,  N  (x  iE 

cooling.  °  1  °>  havc  bccn  neglected  except  in  the  treatment  of  fine-structur 
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Table  3.  Radiative  Processes  Regarded  as  Instantaneous 


Transition 

Initiating  Process* 

r,  sec4" 

rR(S,S0,  eV 

0(3s  V)  -*  0(3P)  +  hi; 

xd  o2(x  3rg+) 

«(-4)* 

9. 14 

0(3s  3S’)  -*  0(3P)  +  hv 

a,  xd  C/2(X  3Eg),  xd  0(3P) 

2. 6(-9) 

9. 52 

0(3p  5P)  -*  0(3s  V)  ♦  hK 

a  02(X  3E  *) 

2. 9(-8) 

1.60 

0(3p  3P)  -*  0(3s  3S°)  +  h v 

a02(X3Eg).  xd  0(3P) 

3.  6(-8) 

1.47 

0+(2s2p4  4P)-  0*(V)  +  hv 

a  0(3P) 

7. 1(-10) 

14.87 

0+(2s2p4  2P)-*  0+(2D3)  +  hK 

.  a  0(3P) 

» 

5. 6(-10) 

23.05 

°2+(A  2nu)  o2+(x  2ng)  +  hv 

a.  1  02(X  3Eg) 

~6.7(-7) 

4.75 

02+(b4Eg’)-.02+(a4nu)  +  hi' 

a,  1  02(X  3Eg) 

<  1.  2(-6) 

2.1 

°a+(*E  ')  -  o2+(A  2nu)  +  h V 

a  02(X  3£g”) 

3.19 

N(3s  4P)  -*  N(2p3  V)  +  hi; 

1 

4 . 3  (-  9) 

10.33 

N(3s  2P)  -  N(2p3  2D*)  +  hi; 

‘  Njt^g") 

1.8(-9) 

8.30 

N(3s  2P)  -*  N(2p3  2P°)  +  hi; 

‘  N2(^g+) 

5. 0(-9) 

7.11 

N(2p4  4P)  -  N(2p3  4S°)  +  hi; 

i  N2(XEg+) 

l-9(-  9) 

10.93 

N(3p  V)  -  N(3s  4P)  +  hi; 

i  N2(XEg+) 

5.  3(-8) 

1.37 

N(3p  4P")  N(3s  4P)  +  hi; 

iNAg+) 

4- 4(-8) 

1.45 

N(3s'2D)  N(2p3  2D*)  +  hi; 

IN/Sg*) 

2.  2(-9) 

9.97 

N(4s  2P)  -  N(2p3  2D°)  +  hi; 

1  *2(\) 

9.  l(-9) 

10.49 

N(3d  2D)  -*  N(2p3  2D*)  +  hi; 

i  N2(1Eg+) 

2.  l(-8) 

10.  65 

Nj(B  3Ilg)  -*  N2(A  3EU+)  +  hi; 

xd  N2(X  XEg+) 

~8.  (-6) 

1.18 

N^a^g)  -  N2(X  ♦  hv 

xd  N2(X  ‘e  *) 

<2.  (-4) 

8.  55 

n2(C  3nu)  N2(B  3ng)  +  hi; 

Xd  N2(X  3Eg+) 

<6.  (-8) 

3.67 

N2(E  3Eg+)  -  N2(A  V)  +  hi; 

Xd  n2(x  *) 

<3.  (-4) 

5. 70 

N2(a'  !£*)  -  N2(X  ^g*)  +  hi; 

xd  N2(X  1Eg+) 

1* 

12.28 

N2(b  -  N2(X  lL+)  +  hi; 

Xd  n2(x  ‘e*) 

da^ 

12.57 

N2<c'  *0  -  N2(X  lj:g+)  + 

xd  N2(X  1Eg+) 

da 

12.93 

N2(b'1Eu+)-N2(X  1Eg*)  +  lw 

xdN2(X  lL*) 

da 

13.21 

N2+(A2nu) -N2+(X  2Eg+)  +  hi; 

a,  1  N2(X  *E  *) 

~1. 2(-5) 

1.0 

N2+(B  2Eu+)-  N2+(X  2Eg+)  +  hi; 

a,  i  N2(X  lL+) 

<7.  (-5) 

3. 17 

*a  h  photoabsorption;  xd  =  collisional  excitation  and  dissociation;  i  =  coilisional  ionization. 

^Monatomic  radiative  iifetimes  from  Wiese  et  al.  (1966);  diatomic  radiative  lifetimes  from  Anderson 
(1971). 

*Deeadic  exponents  are  enclosed  in  parentheses. 

^Estimate;  see  Hcrzberg  (1950). 

^a  =  dipole  allowed,  Chung  and  Lin  (1972). 
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Table  4.  Calculated  Electron  Flux  for  Energy  Groups 


U.,  eV  - 

2  - 1 
Flux,  (cm  sec  eV) 

_ L 

180  km 

200  km 

250  km 

1. 1 

3. 13(9)+ 

2. 11(10) 

7.46(11) 

1.4 

4.26(8) 

2.  85(9) 

1.78(11) 

1.8 

6.38(7) 

3. 17(8) 

2.71(10) 

2.3 

1.07(7) 

5. 10(7) 

2.  09(9) 

2.  9 

1.27(7) 

5.45(7) 

3.  69(8) 

3.6 

4.27(7) 

1.  24(8) 

2.  54(8) 

4.  5 

3.36(7) 

9.  90(7) 

2.  50(8) 

5.7 

2.  52(7) 

7.63(7) 

2.  20(8) 

7.2 

1.46(7) 

4.  90(7) 

1.91(8) 

9.0 

4.  59(6) 

1.  70(7) 

1.03(8) 

11.0 

1.83(6) 

6.83(6) 

5.  53(7) 

14.0 

8.07(5) 

3.02(6) 

2.  66(7) 

18.0 

4.  56(5) 

1.93(6) 

1.  84(7) 

23.0 

3.27(5) 

1.  81(6) 

2.42(7) 

29.0 

2.44(5) 

1.  21(6) 

1. 13(7) 

36.0 

1.  62(5) 

8.  09(5) 

5.  97(6) 

45.0 

1.77(5) 

8. 13(5) 

4.  22(6) 

57.0 

1.23(5) 

4. 17(5) 

1.65(6) 

72.0 

6.61(4) 

1.  26(5) 

'  2.73(5) 

90.0 

4.90(4) 

8.46(4) 

1.  58(5) 

110.0 

3. 17(4) 

4.32(4) 

6. 12(4) 

140.0 

2.02(4) 

2.  63(4) 

3.  54(4) 

180.0 

1.42(4) 

1.74(4) 

2.  23(4) 

230.0 

5.  64(3) 

6.71(3) 

8.  50(3) 

290.0 

1.13(3) 

1.30(3) 

1.61(3) 

360.0 

1. 18(2) 

1.38(2) 

1.75(2) 

Decadic  exponents  are  enclosed  in  parentheses. 


r 


Table  5.  Electron  Production  and  Loss  Rates 


Rate  (electrons  cm~^  sec~*) 
km  200  km  '  250" km 


Photoionization  (PI) 

0(3P) 

■1 

7.6 

29.4 

117.9 

N2<X  X  > 

4.0 

11.8 

38.5 

°2(X  g‘> 

8.3 

6.4 

4.4 

no(x  zn) 

6.8 

3.6 

1.3 

Total 

26.7 

51.2 

162.1 

Collisional  Ionization  (Cl) 

0(3P) 

1 

15. 1 

30.5 

68.2 

N2<X  Sg  > 

13.0 

17.2 

15.9 

°2<X 

1.2 

1.4 

1.0 

Total 

29.  'i 

49.1 

85.1 

Dissociative  Recombination  (DR) 

NO+(X  1L+) 

661.1 

715.4 

894.6 

°2  (X 

N2  <x  sg  > 

203.7 

298.7 

0.9 

224.2 

68.7 

Total 

864.8 

1015.0 

1187. 5 

DR/  (PI  +  Cl) 

15.4 

10.1 

4.8 

t 


t 
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TaMe  6'  25*2™  .Rat?  °f  Formation  of  Species  from  Photoabsorption 
and  Slowing- Down  Photoelectrons  uiudusorpcion 


Relative  Rates,  of  Formation  Relative  Rates  of 

Species  180  to  200  km  250  km  Species  180  ton  200  km  250  km 


0("P) 

0. 14 

0.05 

- 

0(!D) 

0.  26 

0.  23 

0. 16 

O^S) 

0.01 

0.01 

0.  01 

0(3s3S°) 

0.01 

0.02 

0.  03 

0(3s  3P) 

0.02 

0.045 

0.10 

0+(V) 

0.04 

0.08 

0.18 

0+(2D°) 

0.  03 

0.05 

0.11 

0+(2P°) 

0.015 

0.025 

0.06 

°2+(x2ng) 

0.  03 

0.02 

- 

°2+(a  \> 

0.01 

- 

66 


N(4S°) 

0.01 

- 

- 

N2(X  ^Jv 

'  0.18 

0.  20 

0.08 

N2(A3£+) 

0.11 

0. 11 

0.06 

N2(E3ng) 

0.02 

0.02 

0.02 

N2(a  \) 

0.01 

0.01 

0  01 

n2(c  3nu) 

0.01 

0.015 

0.01 

N2+(X  \+) 

0.015 

0.02 

0.035 

N2+(a\) 

0. 03 

0.04 

0.05 

n2+(b  V) 

0.01 

0.01 

0.01 

Table  7.  Relative  Rates  of  Change  in  Energy  Modes  after  Photoabsorption 
and  after  Establishment  of  Electron  Energy  Distribution 


M* 

180  km 

200  km 

250  km 

PA+ 

ssf 

PA 

SS 

PA 

SS 

E 

0.45 

0.04 

0.53 

0.05 

0.  52 

0.12 

• 

K 

C.  06 

0.07 

0.02 

0.03 

-0.0 

-0.0 

• 

R 

0. 10 

0.19 

0. 10 

0.  23 

0.06 

0.  21 

« 

X 

v 

-0.0 

0.01 

0.0 

0. 01 

-0.0 

o 

• 

o 

l 

• 

D 

0. 10 

0.11 

0.03 

0.04 

0.01 

0.01 

• 

I 

0.16 

0.36 

0.22 

0.43 

0.34 

0.  52 

• 

X 

e 

0.13 

0.22 

0. 10 

0.21 

0.07 

0. 14 

*M  =  dM/dt;  E,  electron  translation;  K,  heavy- particle  translation; 

R,  radiation;  Xy,  vibrational  excitation;  D,  dissociation;  I,  ionization; 
xe,  electronic  excitation.  Reference  state  is  N0(X  +)  +  O  (X 

2  g  2  g  '' 


PA  =  photoabsorption;  SS  =  steady  state  in  high-energy  electron  subset. 


Figure  Captions 


1  Differential  electron-energy  spectra  for  electron  collisional  ionization 

of  N2(X  Eg):  Comparison  of  model  in  this  paper  (X  =  0  in  equation  13) 
with  experimental  data  of  Opal  et  al.  (1971). 

2  Electron-flux  spectra  for  a  sample  problem  in  which  the  photoelectrons 
lose  their  energy  by  only  Coulomb  collisions  with  ambient  electrons: 
Comparison  ol  simple  model  with  rigorous  treatment  of  Krinberg  (1973). 
See  Section  2.  2  for  detailed  conditions. 

3  Calculated  and  measured  electron-flux  spectra  at  White  Sands,  N.  M. 
on  February  8,  1971,  for  altitudes  (and  solar  zenith  angles)  of  180  km 
(X  -  90.  2°),  approximately  200  km  (x  ~  90. 1°),  and  250  km  (x  M  90.  0°). 
Smooth  curves  represent  measurements  (Hays  and  Sharp,  1973); 

histogram-like  curves  represent  results  of  calculations.  Note  the 
broken  abscissa. 

4  Calculated  electron-flux  values  in  the  energy  range  from  14  to  36  eV 
at  altitudes  of  200  and  250  km:  Comparison  of  electron-flux  values 
with  (X)  and  without  (O)  the  contributions  from  electrons  formed  by 
absorption  of  the  He  II  303.  8  X  solar  line  by  O(^P)  and  N  (X^E  +) 

2V  g 

The  values  of  the  fluxes  are  plotted  at  the  characteristic  energies  of 
the  energy  groups  (see  Table  1). 
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Calculated  electron-flux  values  In  the  energy  range  from  29  to  90  eV 
at  altitudes  of  180,  200,  and  250  km:  Comparison  of  electron-flux  values 
with  (X)  and  without  (O)  the  contributions  from  electrons  formed  by  photo¬ 
absorption  of  radiation  in  the  wavelength  interval  from  153  to  205  A.  The 

values  of  the  fluxes  are  plotted  at  the  characteristic  energies  of  the  energy 
groups  (see  Table  1). 

Calculated  differential  production  rate  for  photoelectrons  in  the  energy 
range  from  20  to  30  eV  at  200  km  for  *  *  90?  1.  The  ion  products  and 
the  associated  electron  energies  (eV)  for  the  peak  values  are: 

0+(4S°)  -  20.05,  23.33  ,  27.18  ,  30.01;  0+(2D°)  -  20. 02,  23.87,  26.7; 
0+(2P°)  -  22. 18,  25. 01;  N^fX  h+)  -  21. 37,  25.  22,  28.  05; 

N2+(A  2nu)  -  20.37,  24.22,  27.05;  N2+(B  V)  -  22.05,  24  .  88. 
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U'  =  Energy  of  Primary 
DP  =  Degraded  Primary 
S  =  Secondary 
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SECTION  3 


THE  EFFECT  OF  UNCERTAINTIES  IN  RATE 
COEFFICIENTS  ON  DENSITIES  OF  E-REGION  SPECIES 


f 


The  Effect  of  Uncertainties  in  Rate  Coefficients 
on  Densities  of  E-Region  Species 

B.  F.  Myers,  D.  A.  Hamlin, 

M.  R.  Schoonover,  and  J.  I.  Valerio 
Science  Applications,  Inc. ,  La  Jolla,  California  92C37 

A  direct  statistical  method  has  been  applied  to  deter¬ 
mine  the  effect  of  uncertainties  in  rate  coefficients  on  computed 
densities  of  E- region  species.  Uncertainties  in  rate  coeffi¬ 
cients  are  described  in  terms  of  Gaussian  error-law  distribu¬ 
tions  with  a  common  fractional  standard  deviation.  Mean 
densities  are  found  to  be  relatively  insensitive  to  rate- 
coefficient  uncertainties;  fractional  standard  deviations  of 
density  distributions  are,  in  general,  less  than  twice  those 
associated  with  rate- coefficient  distributions,  provided  the 
lati.Tr  fractional  standard  deviations  are  not  greater  than  0.45. 


Many  reactions  are  necessary  to  describe  chemical  processes  in 
the  atmosphere.  In  such  chemical  systems,  it  is  difficult  to  evaluate 
inaccuracies  in  the  computed  species  densities.  These  inaccurares  stem 
from  uncertainties  in  the  input  data  and  from  omissions  or  approximations 
of  relevant  physical  and  chemical  processes.  One  of  the  uncertainties, 
that  associated  with  reaction- rate  coefficients,  is  treated  hereby  a  direct 
statistical  method  for  the  case  of  an  E  region  chemical  reaction  system 


operating  over  short  time  intervals.  The  calculations  presented  here  are 
meant  to  be  illustrative  rather  than  exhaustive. 

Sufficient  information  about  the  errors  and  the  distribution  of  errors 
associated  with  measurements  of  rate  coefficients  is  not  generally  avail¬ 
able.  Accordingly,  in  the  calculations  described  here,  uncertainties  in 
the  rate  coefficients  were  assumed  to  be  describable  in  terms  of  the 

Gaussian  error-law  with  a  mean  value  of  the  rate  coefficient,  k  ,  and  a 

o 

standard  deviation,  s;  the  effects  of  variations  in  s  on  species  densities 
and  their  distributions  were  examined.  For  simplicity,  a  common  frac¬ 
tional  standard  deviation  was  assigned  to  the  rate- coefficient  distributions. 

Consider  a  system  consisting  of  n  reactions  and  corresponding 
rate  coefficients;  let  m(^n)  of  the  rate  coefficients  have  uncertainties 
describable  in  terms  of  the  Gaussian  error-law  distribution.  A  set  of 
rate  coefficients,  jk^},  is  formed  by  randomly  selecting  from  each  k.- 
distribution  for  j  equals  1  to  m.  The  set  of  m  rate  coefficients  and  the 
remaining  n-m  coefficients  are  used  to  calculate  the  species  densities  by 
integrating  the  differential  equations  describing  the  system  over  some 
interval  of  time.  The  selection  of  the  set  jkjj  and  the  integration  of  the 
differential  equations  are  repeated  a  number  of  times  to  obtain  an  indi¬ 
cation  of  the  true  mean  value  and  spread  in  species  densities. 

The  random  values  of  k^  were  sampled  from  the  assumed  distribu¬ 
tion  by  use  of  the  rejection  technique  [Kahn,  1956].  Random  numbers  used 
in  this  technique  were  generated  by  an  adaptation  [McGrath  and  Irving,  1973] 


2 


of  the  standard,  multiplicative  congruential  method  [Coveyou  and 
MacPherson,  1967];  the  adaptation  is  independent  of  the  computer  used 
to  generate  the  numbers.  With  the  use  of  the  Gaussian  error- law  dis¬ 
tribution,  values  of  k.  <  0  occur  in  the  sampling  process;  the  number  of 

J 

such  values  increases  as  s  increases.  These  values  of  k.  were  discarded. 

J 

The  number  discarded  for  the  calculations  described  here  exceeds  0.  05% 

of  the  total  sample  population  only  for  the  largest  value  of  s  used;  in  this 

case,  1.  5%  of  the  sampled  k.  was  discarded.  In  practice,  the  retained 

J 

-2 

values  of  k.  were  equal  to  or  greater  than  10  kQ  (and,  alsoby  choice,  equal  to 
o 

or  less  than  10  k  ). 

o 

The  set  of  reactions  representing  the  E- region  system  is  given  in 
Table  1;  the  system  was  used  under  conditions  applying  to  an  altitude  of 
approximately  140  km.  For  this  system,  n  =  25  and  m  =  17;  the  relative 
values  of  rate  coefficients  for  reactions  1  through  8  and  for  reactions  10 
and  11  were  kept  constant  while  the  absolute  values  of  the  rate  coefficients 
for  reactions  1  and  10  were  varied  in  the  manner  described  above.  The 
differential  equations  for  the  species  densities  were  integrated  over  a 
problem-time  interval  ox  six  minutes  for  each  of  3000  sets  |k.}.  By  six 
minutes  in  problem  time,  initial  transients  were  unimportant.  For  3000 
runs  of  six  minutes  in  problem  time,  2-1/3  minutes  of  CP  time  on  a 
CDC  7600  computer  were  required.  An  additional  calculation  with  the 
set  of  rate  coefficients  listed  in  Table  1  was  also  made.  The  initial 
species  densities  are  listed  in  Table  2  for  the  starting  time  of  the 


3 


calculation,  1200  hrs.  For  each  set  of  3000  calculations,  the  error-law 
distribution  of  rate- coefficient  values  was  characterized  by  the  mean 
value  of  the  rate  coefficient,  kQ,  referred  to  below  as  the  nominal  value, 
and  by  a  fractional  standard  deviation,  fgk  =  s/kQ,  of  0. 15,  0. 30,  and 
0.45  in  successive  sets  of  calculations. 

The  results  of  the  calculations  are  summarized  in  Table  3.  Consider 

first  the  group  of  species  consisting  of  the  charged  particles  plus  N(2D°). 

The  following  features  are  noted  for  this  group:  (1)  The  average  species 

densities  depart  by  less  than  10%  from  the  species  densities  calculated  by 

using  the  set  of  nominal  rate  coefficients,  provided  f  so.  30.  (2)  The 

fractional  standard  deviations  for  species  densities,  f  ,  are,  roughly, 

between  one  and  two  times  f^.  Furthermore,  fsn  increases  monotonically 

with  decrease  in  density,  n,  except  for  N_+.  The  increase  in  f  is 

*  sn 

roughly  a  factor  of  two  for  a  decrease  in  n  from  about  105cm"3to  about  101  cm"3. 
(3)  The  density  distributions  have  pronounced  positive  coefficients-of- 
skewness,  f},  which  increase  with  increasing  f  . 

The  source  of  the  positive  skewness  can  be  appreciated  with  a 
simple  example.  Consider  the  species  densities  to  be  in  a  steady-state 
condition,  a  reasonable  assumption  for  the  present  calculations.  Let  the 
nominal  rate  coefficients  associated  with  production  processes  for  a  given 
species  be  multiplied  by  the  factor  (1  +  6)  while  those  associated  with  the 
loss  processes  of  the  given  species  be  multiplied  by  the  factor  (1  -  €). 

Then  the  nominal  value  of  the  density  of  the  given  species,  n  ,  will  be 


multiplied  by  the  factor  F  =  (1  +  e)/(l  -  c).  In  the  inverse  case,  the 
nominal  value  of  the  density  will  be  multiplied  by  F  1.  On  a  linear  scale 
of  density,  n  F~*  will  be  closer  to  n  than  n  F  is  to  n  :  hence,  the 

J  O  0  0  O 

positive  skewness  obtains. 

The  increases  in  f  and  0  with  increasing  f  .  ,  as  shown  in 

sn  sk 

2 

Table  3  for  the  group  consisting  of  charged  particles  plus  N(  D°),  is 
illustrated  in  Fig.  1  where  a  typical  set  of  histograms  is  drawn;  the  species 
represented  is  N0+.  The  class  intervals  in  these  differential  frequency 
histograms  are  equal  to  0.  2  of  the  standard  deviation  of  the  density  dis¬ 
tribution  for  each  of  the  three  f  k  cases.  (Note  that  the  means,  standard 
deviations,  and  coefficients  of  skewness  in  Table  3  were  calculated  from 
the  original  set  of  3000  densities  prior  to  constructing  the  histograms. ) 

The  mean  density  value  for  each  distribution  is  indicated  by  the  placement 
of  the  associated  arrow.  One  sees  that,  as  f  k  increases,  the  mean  density 
and  the  spread  of  the  distribution  increase.  Also,  the  large  values  of  0 
evidently  arise  from  a  small  number  of  calculations  yielding  high  densities. 
For  the  group  of  species  consisting  of  charge  particles  plus  N(  DJ),  such 
calculations  which  yield  densities  beyond  2.  5  standard  deviations  from  the 
mean  density  represent  less  than  3%  of  the  total  number  of  calculations  for 
f  =  0.45;  for  the  smaller  values  of  f  , ,  such  percentages  are  much 
smaller.  When  calculations  of  the  standard  deviation  and  coefficient  of 
skewness  were  made  with  densities  lying  within  2.  5  standard  deviations  of 
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the  previously  computed  means,  then  for  fgk  =  0.  30,  the  values  of  s  were 
reduced  by  only  12  to  24%  whereas  values  of  ,8  were  decreased  by  65  to  90%. 

For  the  group  consisting  of  the  last  five  species  in  Table  3,  the 
following  features  occur:  (1)  The  average  densities  are  he  same  as  those 
calculated  by  using  the  set  of  nominal  rate  coefficients.  (2)  The  fractional 
standard  deviations  for  species  densities  are  small  but  increase  with 
mcreasmg  fgk.  (3)  The  skewness  of  the  density  distributions  is  small. 

The  above  calculations  were  repeated,  for  f  =  0. 30,  by  using  a 
different  sequence  of  random  numbers  from  the  random  number  generator. 
Use  of  the  new  sequence  did  not  result  in  any  essential  changes.  To  obtain 

an  indication  of  the  effect  of  restricting  the  k  to  values  equal  to  or  greater 

-  2  ^ 
than  10  kQ,  the  above  calculations  were  repeated  by  using  the  subset  of 

reactions  4,  19,  and  20  which  is  not  coupled  to  the  reaction  system  (as¬ 
suming  constant  densities  for  O,  <>2,  and  Ng)  but  yet  describes  completely 
the  processes  involving  one  of  the  species,  0+(2D°).  The  changes  occur¬ 
ring  with  limits  of  k.  given  by  kQ(l  ±  0.  99),  rather  than  kQ(l  +  99)  and 

k0(l  "  °*  99)  as  above’  were  not  larger  than  in  the  cases  of  using  different 
sequences  of  random  numbers, 

For  simplicity  in  these  calculations,  the  problem  of  assessing  the 
errors  given  to  experimental  values  of  rate  coefficients  has  been  circum¬ 
vented.  If  one  consults  a  listing  of  such  errors  for  the  chemical  reactions 
of  Table  1,  as  for  example,  Bortner  and  Baurer  [1973],  and  assumes  that 
the  errors  given  represent  one  standard  deviation  in  a  Gaussian  distribution 
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of  rate- coefficient  values,  then,  for  the  current  state  of  knowledge  about 
the  accuracy  of  rate  coefficients,  the  best  representative  for  a  calculation 
of  the  above  type  would  be  one  using  0.  30  <  f  <  0.  45.  Note  that  if  the 
assumption  were  acceptable,  a  more  reliable  calculation  could  be  mad,' 
with  the  present  method  since  different  fgk  can  readily  be  assigned  to 
different  rate  coefficients. 

The  direct  statistical  method  used  above  has  the  disadvantage  that 
for  long  problem-times,  the  method  is  expensive.  Yet,  the  distributions 
of  species  densities  will  change  with  time  so  that  the  time  dependence  can 
be  important.  The  above  application  is  also  limited  by  use  of  the  Gaussian 
error-law  distribution  which  will  become  a  poorer  representation  as  the 
standard  deviation  in  the  rate- coefficient  distribution  increases;  this 
limitation  can  be  readily  overcome  by  using  a  more  appropriate  distribu¬ 
tion  function.  Finally,  the  above  method  does  not  readily  give  information 
about  the  relative  importance  of  rate  coefficients.  For  the  interested 
reader,  reference  is  given  to  the  related  work  of  Cukier  et  al.  [1973]  and 
Schaibly  and  Shuler  [1973]  in  which  a  method  is  given  for  identifying  those 
rate  coefficients  which  significantly  affect  species  densities. 
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Tabic  1 


Processes  and  Nominal  Values  for  Rate  Coefficients 


Reaction 

No. 

Processa 

4o> 

Reference 

1 

Ng  +  ht»  -*  Ng+  +  e 

3. 0(-8) 

c 

2 

Og  +  hy  -*  Og+  +  e 

7. 0(-8) 

c 

3 

O  +  hv  O4"  +  e 

1.  2(-8) 

c 

4 

O  +  hi'  -»  0+(2D° )  +  e 

1.3  (-8) 

c 

5 

NO  +  hy  -»  NO+  +  e 

7.  8(-7) 

c 

6 

Ng  +  hv  •*  2N 

4.  6(-10) 

c 

7 

Og  +  hi/  ->  20 

1.7(-7) 

c 

8 

NO  +  hv  -»  N  +  O 

2. 3(-7) 

c 

9 

02+  +  e  -*  0  +  O 

1.2(-7) 

d,  e 

10 

NO+  f  e  -*  N  +  O 

4. 4(-8) 

U  g 

11 

NO+  +  e-*  N(2D°)  +  0 

1.8(-7) 

U  g 

12 

N  +  +  e •*  N(2D°)  +  N 

1- 3(-7) 

e,  h 

13 

n2+  +  o2-n2  +  o2+ 

3.  O(-ll) 

i>  j 

14 

N?+  +  0  -*  NO4  +  N 

1. 4  (- 10) 

k,  1 

15 

o+  +  o2  -» o  +  o2+ 

1.  5(- 11) 

i,  m 

16 

0+  +  N2->  NO+  +  N 

6. 3(- 13) 

i,  m,  n,  o 

17 

02+  +  NO->  Og  +  NO+ 

6.  3(-10) 

j>  k 

18 

02+  +  N  -*  NO+  +  0 

1.  8(- 10) 

P 

19 

o+(2d°)  +  o2^o  +  o2+ 

3.  6(- 10) 

q 

20 

o+(2d°)  +  n2-*o  +  n2+ 

4.  5(-10) 

r 
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Table  1  (Continued) 


Reaction  for 

No- _ Process3- _ 1200  hrb  Reference 

21  N  +  Og  -»  NO  +  O 

22  N(2D°)  +  C>2  -*  NO  +  O 

23  N(2D°)  +  O-*  N  +  O 

24  N  +  NO  -*  N2  +  O 

25  N(2D°)  +  NO-*  N2  +  O 

a  _  '  ‘  “  - - - - 

Only  ground  electronic  states  are  involved  unless  indicated  otherwise. 

The  neutral  particle  temperature  was  taken  from  the  CIRA-65  Model- 5 
atmosphere  to  be  610'K  for  1200  hrs;  units  of  ko  are  sec"1  for 
processes  1  to  8  and  cm3  sec-1  for  the  remaining  processes.  The 
decadic  exponent  is  enclosed  by  parentheses. 

Calculated  with  solar  flux  routine  using  flux  values  and  photon  group 
structure  recommended  by  Hinteregger  [1970],  for  conditions  of  1200  hrs 
at  140  km  with  solar  zenith  angle  of  21.  6°. 

^Cunningham  and  Hobson  [1972]. 

eMehr  and  Biondi  [1969]. 

^Weller  and  Biondi  [1968]. 

^Walls  and  Dunn  [1973]. 

hMichels  [1972]. 

1Dunkin  et  al.  [1968]. 

JJohnsen  et  al.  [1970]. 

kFehsenfeld  et  al.  [1970]. 

^eynaber  et  al.  [1973]. 

mJohnsen  and  Biondi  [1973]. 

nBreig  et  al.  [1973]. 

°Van  Zandt  and  O'Malley  [1973]. 

PGoldan  et  al.  [1966]. 

A.  Rutherford  (unpublished  data,  1973). 


4.  3(- 14) 
0.  9(- 12) 
1.  0(-12) 
2.  2(- 11) 
1.  3(-10) 


s 

t,  u 

no  data; 
estimated 

v 

t,  w 
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r 


r  Rutherford  and  Vroom  [1971]. 
sClark  and  Wayne  [1970]. 

*Lin  and  Kaufman  [1971 J. 
uSlanger  et  al.  [1971]. 
vSchif:  [1969]. 
wBlack  et  al.  [1969]. 
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Table  2 


Initial  Densities  of  Species  at  t  =  1200  hrs* 


Species 

Density, 

cm"^ 

Species 

Density, 

cm'^ 

N2 

6.4(10) 

N2+ 

6.7(2) 

°2 

9.7(9) 

°2 

2.6(4) 

0 

2.2(10) 

0+ 

1.  2(3) 

NO+ 

5.0(7) 

0+(2d°) 

8. 0(0) 

N 

1.  9(6) 

NO+ 

8. 1(4) 

N(2D) 

3.3(4) 

* 

e 

« 1.1(5) 

t  =  starting  time  of  calculation, 
t  ° 

Densities  from  Strobel  [1972]. 

In  the  calculations,  the  electron  density  was  set  equal  to 
the  exact  sum  of  ion  densities. 


Table  3 

Summary  of  Results  on  Propagation  of  Rate-Coefficient  Uncertainties  for  E-Region  Densities 

fsk  =  °-15  fsk  =  0.  30  fsk  =  0.45 

Specie  nQ  %/<">  I  fsn  ?  no/<n>  I  fsn  ^  no/<n>  I  fSn  I  5 


density  of  species  computed  with  nominal  rate-coefficient  set. 

average  density  of  species  computed  with  3000  sets  of  randomly-selected  rate  coefficients. 

sx/<x>,  where  sx  is  the  standard  deviation  of  the  x-distribution  and  x=k  for  rate- coefficient 
distribution  and  x=n  for  density  distribution  of  species. 

skewness  coefficient  =  third  central  moment/ (standard  deviation)^ 


I 


1  ig.  1.  Differential  frequency  histograms  for  NO+  in  cases  of  f  =  o.  15 
0.30,  and  0.45.  The  class  intervals  are  equal  to  0.  2  times  the  standard 
deviation  for  each  of  these  cases.  Note  the  change  of  scales  for  the  lower 
panel  in  which  the  histograms  are  extended  above  densities  of  2  x  lo5  cm'3 

The  mean  density  value  for  each  distribution  is  indicated  by  the  placement 
of  the  associated  arrow. 
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SECTION  4 


IONOSPHERIC- ELECTRON  HEATING 
BY  DISSOCIATIVE  RECOMBINATION 


Research  Note 


Ionospheric-Electron  Heating  by  Dissociative  Recombination 

Q  .  A  Dl  and  B.  F.  Mvers 

Scence  Apphcations,  Inc. ,  La  Jolla,  California  92037,  U.  S.  A. 


Abstract  -  It  is  show,  that  the  loss  of  ionospheric  electrons  by 
dissociative  recombination  effectively  heats  the  electron  gas  at 
a  rate  which  is  a  few  per  cent  of  the  main  heating  rate  due  to 
Photoelectrons  but  which  exceeds  some  of  the  rates  for  minor 
cooling  processes  frequently  included  in  energy-balance  studies. 


A  number  of  effects  have  been  considered  in  the  literature  for 
heating  and  cool.ng  the  ionospheric  electron  gas  (Stubbe  and  Vamum,  1972). 
In  this  Note  we  assess  an  overlooked  effect  of  intermediate  importance. 

In  the  E  and  F  regions,  electrons  are  lost  mainly  by  dissociative 
recombination  with  NO+  and  Og+.  When  an  electron  is  lost,  its  kinetic 
energy,  E,  is  removed  from  the  thermal  energy  of  the  electron  gas  at 
perature  8(eV).  The  electron  gas  is  effectively  heated  if  E  <  39/2 
and  cooled  if  E  >  39/2.  The  effective  energy-return  to  the  electron  gas 
per  recombination  with  a  tfven  molecular  ion,  e,  is 


«  =  1  0  -  <E>  . 


(1) 


1 


Here,  for  a  dissociatively- recombining  electron  of  velocity  v  ~^2E/m 
and  cross-section  o(E),  its  mean  energy,  (E>,  is  given  by 


E  o(E)  v  f(E) 


dE  , 


(2) 


where  f(E)  is  the  Maxwellian  energy-distribution  function,  and  the  rate 
coefficient  corresponding  to  a(E)  is 


k  h  <cv>  =  f  a(E)  v  f(E)  dE  .  (3) 

J0 

The  heating  rate  for  the  electron  gas  is 
t 

*  dH  o  i 

H  H  *df  =  k  [e]  [i]  e  ,  eV  cnn  sec"1  (4) 

where  [e]  is  the  electron  density  and  [i]  is  the  molecular-ion  density. 

The  dependence  of  c  on  E  or  of  k  on  6  is  not  well  established 
for  ionospheric  conditions.  Thus,  we  shall  consider  a  range  of  possibilities 
For  simplicity,  we  assume  that  the  cross  section  can  be  expressed  as  a 
power  law, 

a  =  oQ  E"n  ,  n  <  2  .  (5) 

Evaluation  of  Equations  (2)  and  (3)  with  use  of  Equation  (5)  leads  to  the 
results 


<E>  _  r(3  -  n)  - 

0  r(2  -  n)  2  ”  n  »  r\<  2  (6) 


and 


k  =  2oQy]2/vm  r(2  -  n)  0-p  , 


(7) 


2 


where  r(x)  is  the  gamma  function  and 


Use  of  Equations  (1),  (6),  and  (8)  gives 

c  =  p0  .  (9) 

Evaluation  of  Equations  (6),  (8),  and  (9)  for  several  values  of  n  gives  the 
results  in  Table  1.  The  experimental  and  theoretical  indications  for  NO+ 
and  02+  are  that  ~  s  p  s  1  and  hence  that  s  c/e  ^  1. 

As  a  practical  matter,  however,  the  dissociative- recombination 
heating  rate  may  be  insensitive  to  the  particular  power- law  rate  coeffi¬ 
cient  used.  For  example,  consider  rate  coefficients  kj  =  kQ1  0  ?1  and 
“P2 

^2  =  ^02  ®  »  define  0^  and  0^  to  be  the  temperatures  for  which  the 

respective  equalities  k^  =  k^  and  k^pj  =  kgP2  are  sa^*s^ed>  and  define  R 

to  be  the  ratio  of  the  heating  rates,  R  =  Hj/Hg  =  kjPj/kgPg.  If  px  <  p2, 

then  R  <  1  for  0  <  0.  but  R  increases  for  9  >  0.  until  R  =  1  at  0  =  0 

K  k  kp 

and  R  >  1  for  0  >  0^  .  Thus,  if  the  problem  of  interest  is  one  for  which 
the  temperature  has  a  small  range  about  0^,  then  the  heating  due  to  dis¬ 
sociative  recombination  will  depend  little  on  the  choice  made  for  k. 

In  Fig.  1  we  illustrate  the  magnitude  of  dissociative-recombination 
heating  relative  to  other  heating  and  cooling  mechanisms  for  ionospheric 
electrons.  The  usual  heating-  and  cooling- rate  curves  in  Fig.  1  are  those 
computed  by  Timothy  et  al.  (1972)  for  their  rocket  flight.  To  compute  the 


3 


curve  for  dissociative- recombination  heating,  we  used  the  densities  and 
temperatures  given  by  Timothy  et  al.  (ib7?):  values  of  [e],  [NO+],  and 
[Og*]  were  taken  from  Table  11  (temperature-condition  c);  values  of  0 
were  taken  from  Fig.  37  (measured  values  for  h>  150  km  and  calculated 
values  for  h  150  km).  We  used 


k(NO+)  =  4.  0  x  l(f7  (T/300)-0,  8 

(10) 

k(02+)  =  2. 1  x  10"7  (T/300)”0, 7 

(ID 

for  the  dissociative- recombination  rate  coefficients  (Bortner  et  al. ,  1973) 
and  Equation  (4)  for  each  of  the  heating- rate  contributions  from  NO+  and 
0^‘  We  see  that  the  heating  rate  due  to  dissociative  recombination  is  a 
small  fraction  (from  2.  6  to  5. 3  per  cent)  of  the  main  heating  rate  due  to 
photoelectrons  but  that  it  does  ex:eed  some  of  the  rates  for  minor  cooling 
processes  frequently  included  in  energy-balance  studies. 

After  completing  this  work,  we  noted  a  comment  by  Dalgarno 
(1973)  that  relates  closely  to  the  topic  discussed  here,  though  it  was  made 
in  an  astrophysical  context  with  respect  to  radiative  recombination:  ’’The 
cold  electrons  recombine  more  readily  and  recombination  is  itself  a 
mechanism  for  increasing  Tg  [the  electron  temperature].  ” 

Acknowledgement  —  This  work  was  partially  supported  by  the  Defense 
Advanced  Research  Projects  Agency  of  the  Department  of  Defense  and  was 
monitored  by  ONR  under  Contract  No.  N00014-72-C-0292. 
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TABLE  1.  Parameters  for  Effective  Energy-Return 
to  Electron  Gas  for  Power-Law 
Dissociative-Recombination  Cross- 
Section 


n* 

p+=  c/e 

<E>/e 

1 

"2 

-l 

5 

2 

0 

l 

2 

1 

3 

0 

3 

2 

1 

l 

1 

3 

3 

l 

1 

2 

^Exponent  in  Equation  (5). 


^Exponent  in  Equation  (7). 
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ELECTRON  COOLING  AND  HEATING  RATES  (eV  cm'3  sec'1) 


Fig.  1.  Calculated  cooling  and  heating  rates  for  ambient  electrons. 
Cooling  rates  (Timothy  et  ab  ,  1972;  Fig.  38):  (1)  elastic 
(positive  ions);  (2)  elastic  (neutrals);  (3)  rotational 
excitation;  (4)  fine-structure  excitation;  (5)  electronic 
excitation;  (G)  vibrational  excitation.  Heating  rates: 

O  photoelectrons  (Timothy  etab  ,  1972,  Fig.  27); 

□  dissociative  recombination  (see  text) 
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SECTION  5 

t.t^J^MPEFlAT^RE  DEpENDENCE  FOR  DISSOCIATIVE 
RECOMBINATION  OF  NO+  IN  E-  AND  F-Rij  JION  MODELS 


Temperature  Dependence  for  Dissociative  Recombination 
of  NO+  in  E-  and  F-Region  Models^ 


B.  F.  Myers 

Science  Applications,  Inc. 

La  Jolla,  California  92037,  U.  S.  A. 


Abstract  -  A  discussion  is  given  which  indicates  that  the  temperature 
dependence  for  the  dissociative  recombination  of  NO+  customarily  used 
in  E-  and  F-region  models  is  not  correct. 


Recent  atmospheric  models,  such  as  Torr  and  Torr  (1969), 
Keneshea  et  al.  (1970),  Stubbe  (1970),  Nisbet  (1971)  and  Riister  (1971), 
which  include  a  description  of  chemistry  in  the  E-  and  F-regions,  use  a 
temperature*  dependence  of  T  or  T  ^  for  the  rate  coefficient  of  dis¬ 
sociative  recombination  of  NO+.  On  the  basis  of  present  experimental 
and  theoretical  evidence  (Bardsley,  1970;  Bardsley  and  Biondi,  1970), 
these  temperature  dependencies  are  incorrect  for  temperatures  to  be 
encountered  in  the  E-  and  F-regions. 

The  appropriate  temperature  dependence  of  the  rate  coffficient 
should  be  that  for  the  dissociative  recombination  of  NO^  which  is  essen¬ 
tially  in  the  ground  vibrational  state.  Although  the  temperatures  in  the 
E-  and  F-regions  correspond  only  to  a  small  degree  of  vibrational  exci¬ 
tation  of  NO  under  the  assumption  of  equilibration  of  energy  modes,  this 
is  not  a  sufficient  basis  for  the  preceding  recommendation.  The  forma¬ 
tion  of  NO+  is  principally  by  way  of  the  processes 


The  temperature  is  specifically  the  electron  temperature  under  condi¬ 
tions  of  nonequilibrium  among  the  translational  modes  of  electrons,  ions 
and  neutral  particles. 

^Published  in  Journal  of  Atmospheric  and  Terrestrial  Physics,  1973, 

Vol.  35,  pp.  1903- 1 904.  Pergamon  Press. 
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(1) 


N2+  +  O  -  N0+  +  N 

02+  +  NO  -  02  +  NO+  (2) 

0+  +N2-  NO+  +  N  .  (3) 

Reactions  (1)  and  (2)  have  large  rate  coefficients  and  are  highly  exothermic. 
Consequently,  the  product  NO+  may  be  formed  with  energy  in  the  vibra¬ 
tional  mode  much  in  excess  of  that  corresponding  to  E-  and  F-region 
temperatures.  This  possibility  cannot  be  dismissed  on  the  basis  of 
present  experimental  results  and  for  reaction  (2)  current  knowledge  about 
charge  transfer  processes  involving  molecular  species  supports  the  possi¬ 
bility  (Rutherford  et  al.  ,  1971;  Tomcho  and  Haugh,  1972).  For  reaction  (3), 
which  has  a  smaller  rate  coefficient  and  is  less  exothermic  than  either 
reactions  (1)  or  (2),  recent  experiments  (Neynaber  and  Magnuson,  1973) 
indicate  that  for  most  of  the  NO+  produced,  75  per  cent  or  more  of  the 
available  energy  (reaction  exothermicity,  vibrational  energy  of  N2  and  rela¬ 
tive  kinetic  energy  of  reactants)  may  be  transferred  to  internal  modes  of 
NO  .  If  NO+  is  formed  with  substantial  energy  in  vibrational  modes,  and 
if  there  were  no  way  to  remove  this  energy  before  dissociative  recombina¬ 
tion,  the  rate  coefficient  to  be  used  could  be  significantly  different 
(O'Malley,  1969;  Bardsley,  1970)  from  the  case  with  energy  in  the  vibra¬ 
tional  modes  corresponding  to  E-  and  F-region  temperatures.  Thus,  in 
the  absence  of  definitive  experimental  results  for  the  distribution  of 
energy  among  the  products  of  reactants  (1),  (2)  and  (3),  and  of  redistribution 
of  energy  involving  vibrational  modes  of  NO+  prior  to  dissociative  recom¬ 
bination,  no  recommendation  could  apparently  be  made  concerning  the 
temperature  dependence.  However,  in  spite  of  the  possibility  for  forming 
NO  in  vibrationally  excited  levels,  when  undergoing  dissociative  recom¬ 
bination  the  NO+  will  be  in  the  ground  vibrational  state,  on  the  average, 
as  a  result  of  radiative  decay.  Radiative  decay  of  vibrationally  excited 
NO+  molecules  occurs  in  the  infrared  spectral  region  with  frequencies  of 


2 


the  order  of  102  sec"1  (Herzberg,  1950;  Miescher,  1956;  Penner,  1959; 
Stair  and  Gauvin,  1967);  consequently  radiative  decay  will  predominate, 
even  ignoring  energetic  and  transition- probability  considerations,  at  ^ 
altitudes  above  about  120  km,  where  the  collision  frequencies  (<  10  sec" 
are  smaller  than  radiative  decay  frequencies. 

Just  what  temperature  dependence  should  be  used  in  E-  and 
F- region  models  remains  to  be  firmly  established.  Perhaps  a  T 
dependence  indicated  by  theory  (Bardsley,  1970)  could  be  used  ir.  the 
interim. 
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